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FOREWORD 


This  report  was  prepared  by  the  McDonnell  Aircraft 
Company,  (MCAIR) ,  a  Division  of  the  McDonnell  Douglas  Corpora¬ 
tion,  St.  Louis,  Missouri  for  the  Aerospace  Dynamics  Branch, 
Vehicle  Dynamics  Division,  Air  Force  Flight  Dynamics  Laboratory, 
Wright-Pat  terson  Air  Force  Base,  Ohio  under  Contract.  AF 
33615-70-C-1272  ,  Project  No.  1370  ,  "Dynamic  Problems  in 
Military  Flight  Vehicles",  Task  No.  137008,  "Prediction  of 
Dynamic  Loads  in  Military  Aircraft."  Mr.  Verner  J.  Johnson 
was  the  Project  Engineer  at  the  Air  Force  Flight  Dynamics 
Laboratory. 

This  program  represents  efforts  by  MCAIR  personnel  of 
the  Engineering  Technology  and  Laboratory  Flight  Test  Divisions, 
as  well  as  by  Dr.  John  Houblot,  Consultant.  The  major  portion 
of  the  technical  effort  was  provided  by  Mr.  Robert  E.  Mullans, 
Senior  Structural  Dynamics  Engineer,  under  the  direction  of 
Dr.  Clark  Lemley,  Project  Dynamics  Engineer.  Dr.  Norman 
Zimmerman  provided  overall  program  management.  The  authors 
wish  to  acknowledge  the  participation  of  Dr.  Houbolt  who  made 
significant  contributions  to  Sections  5  and  6  of  this  Final 
Report.  The  manuscript  was  released  by  the  authors  for  publi¬ 
cation  in  March  1972. 

This  technical  report  has  been  reviewed  and  is  approved. 


WALTERS  J.  MYKYTOW 

Asst,  for  Research  and  Technology 
Vehicle  Dynamics  Division 
AF  Flight  Dynamics  Laboratory 


ABSTRACT 

This  report  describes  an  empirical/theoretical  study  program  to  develop 
a  first  generation  method  for  predicting  the  buffet  dynamic  loads  and  fighter 
aircraft  responses  that  occur  during  transonic  maneuvers.  The  prediction 
method  applies  to  wing  upper  surface  loads  acting  on  the  airplane  and  the 
induced  aircraft  dynamic  responses.  Wind  tunnel  tests  of  a  10#  scale  model 
of  the  YF-1*E  aircraft  were  used  to  enhance  understanding  of  the  buffet  mech¬ 
anism  for  lifting  surfaces  and  to  obtain  data  needed  to  deduce  the  surface 
load  formulas.  A  first  test  phase  (Series  l)  gave  oil  flow  and  tuft  infor¬ 
mation  regarding  flow  on  the  upper  wing  surfaces.  A  second  test  phase  (Series, 
II)  was  used  to  acquire  both  steady  mid  fluctuating  pressure  and  strain  data 
as  well  as  acceleration  response. 

The  fluctuating  pressure  data  were  converted  to  power  and  cross-power 
spectra  and  a  mathematical,  model  was  formulated  to  describe  the  amplitude 
and  the  spatial  and  spectral  distributions  of  the  pressures  on  the  wing. 

The  mathematical  pressure  distribution  was  then  used  as  the  forcing 
function  in  a  dynamic  response  computer  program  to  assess  the  structural 
response.  Model  pressure  data  and  the  predicted  airplane  responses  were  in 
reasonable  agreement  with  flight  test  data, 

Tn  the  process  of  developing  the  method  for  buffet  loads  prediction,  a 
great  deal  of  valuable  information  was  collected  pertaining  to  the  nature 
of  the  buffeting  flow  over  a  lifting  surfaee  mid  the,  resulting  (structural 
responses,  'litis  data  together  with  Its  interpretation  are  included  in  this 
report  mid  supplementary  iv|H»rttt  that  aiv  referenced  herein. 
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for  both  positive  and  negative  frequencies.) 

Nondimensional  distance 

One-half  period  of  time  over  which  spectral  functions  are 
determined,  theoretically  approaching  infinity 


Free  stream  velocity 

Convection  velocity 

Distance  aft  of  wing  leading  edge 

Chordwise  wing  coordinate 

Distance  outboard  from  aircraft  centerline 

Spanwise  wing  coordinate 

Reciprocal  of  complex  frequency  response  function 
Displacement 

Streamwi se  decay  coefficient 
Wing  angle  of  attack  at  buffet  onset 
Wing  angle  of  attack 
Spanwise  decay  coefficient 

Magnitude  of  cross  power  spectral  density  function  normalised  by 
dividing  by  square. root  of  product  of  two  corresponding  power 
spocJ  4,1  density  functions,  Also,  square  root  of  coherence 
function.--  -  . 

Spatial  decay  coefficient 

Ratio  of  damping  to  critical  damping 

Ratio  of  the*  peak  level  to  the  low  frequency  level  in  a  power 
spectral  density  function 

Phase  angle  of  cross  power  spectral  density  function 
Parameter  in  mathematical  model  for  power  spectral  density 
Streamwiae  distance 
Spanwise  distance 

Autocorrelation  function  of  u  pressure  nonnaliaed  by  dividing 
by  the  mean  square  value  of  the  pressure 
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P  ( i' * t  )  Cross-correlation  function  of  two  pressures  normalized  by- 

dividing  by  the  product  of  the  rms  values  of  the  pressures 

a  Root  mean  square  value  of  a  quantity  indicated  by  a  subscript 

or  the  context 

x  Time  shift  for  autocorrelation  or  cross-correlation  functions 

(j>  Modal  deflection 

ip  Power  spectral  density  function  defined  only  for  posxtive 

frequencies  and  normalized  by  dividing  by  the  corresponding 
mean  square  value 

a)  Frequency  in  radians  per  second 

id  Frequency  parameter  in  mathematical  model  for  pressure  auto¬ 

correlation  function 

Subscripts  and  Superscripts 

a,  b  Indicate  geometrical  points  on  structure 

d  Indicates  modified  by  decay  correction  factor 

F  Indicates  filtered  time  function 

FS  Indicates  full  scale 

i,  j  Indicate  modes  i  and  j,  respectively 

k,  &  Indicate  areas  k  and  ,  respectively 

M  Indicates  relationship  to  bending  moment 

M  Indicates  model  {as  opposed  to  full-scale) 

p  Indicates  relationship  to  pressure 

q  Indicates  relationship  to  normal  coordinate 

z  Indicates  relationship  to  displacement 

l,  g  Indicate  relationship  to  arbitrary  functions  1  and  2 

Matrix  Notation 


Square  matrix 
Row  matrix 

Transpose  (of  row  matrix,  for  example) 
Column  matrix 


Ml  sceX-Uuteous  Notation 

—  Denotes  mean  v&luu 

*  Denotes  complex  conjugate 

r/  ff  ■uentlation  of  a  variable  with  respect  to  time  is  indicated  by  dots 
over  the  variable 
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INTRODUCTION 

Airplane  wing  buffet  occurs  when  an  aircraft  is  in  a  high  angle  of 
attack  flight  regime  and  the  air  flow  on  the  upper  wing  surface  becomes 
detached.  Air  flow  in  the  detached  regions  becomes  turbulent  thereby  giving 
rise  to  fluctuating  pressures  on  the  wing,  and  the  possibility  of  induced 
pressures  on  the  tail.  While  some  loss  of  lift  may  be  experienced,  there  is 
an  appreciable  range  of  angles  of  attack  beyond  the  buffet  boundary  in  which 
an  airplane  can  maneuver.  The  upper  limit  is  usually  determined  by  wing  rock 
or  stall.  This  report  describes  the  development  of  a  first  generation  buffet 
load  prediction  method  for  a  modern  fighter  aircraft  in  the  transonic  flight 
regime . 

1.1  History 

Buffet  investigations  that  have  dealt  with  aircraft  can  be  divided  into 
two  distinct  time  periods  -  those  that  occurred  prior  to  the  intensive  U,S. 
space  effort  (which  dominated  the  late  1950's  and  196Q's),  and  those  which 
have  occurred  since  about  I96B. 

1.1.1  Earlier  Studies  -  Impetus  for  the  earlier  studies  was  provided 
mainly  by  the  necessity  to  understand  the  nature  and  severity  of  the  buffeting 
loads  and  their  potentially  destructive  effect  on  aircraft  structures .  While 
the  results  of  buffet  eould  be  readily  noted  by  ride  quality  and  sometimes 
by  aircraft  damage  (Reference  l) ,  analytical  attack  on  the  problem  was  ham¬ 
pered  by  the  transient  and  random  nature  of  the  airloads  causing  buffet. 

The  buffet  problem  eunnot  be  bundled  in  a  deterministic  manner.  While  mathe¬ 
matical  theory  was  available,  it  remained  for  the  speed  and  capacity  of  modern 
electronic  computers  and  signal  processing  equipment  to  open  the  door  to  a 
more  thorough  treatment  of  the  theoretical  aspects  of  the  problem. 

1.1. ‘i  Recent  Studies  -  Recent  studies  have,  in  a  sense,  been  a  continua¬ 
tion  of  the  earlier  efforts.  However,  the  scope  of  needed  buffet  investiga¬ 
tions  has  been  greatly  increased  by  the  desire  to  utilise  aircraft  maneuver 
capabilities  beyond  buffet  onset  even  well  into  the  transonic  flight  regime. 

The  Air  Force  (References  U ,  3,  and  k)  and  the  NAHA  (Reference.,  5  and  6) 

.xv it  conducted  in-depth  investigations  of  buffet  both  analytically  and  oxporimeu 
tally  and  Reference  7  indicates  some  of  the  effort  being  carried  on  in  Great 
rr 1  tain.  Most  of  these  studies,  however,  have  been  concerned  mainly  with 
definition  of  buffet  unset ,  and  with  the  measurement  of  static  pressures  and 
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loads.  The  use  of  wide  band  fluctuating  pressure  sensors  to  measure  aircraft 
buffet  pressures  is  still  very  limited,  and  most  effort  in  this  direction  has 
been  sponsored  by  NASA  Ames  and  by  the  Air  Force  (as  described  in  this  con¬ 
tract  ) . 

Although  some  of  the  technology  gained  from  space  launch  vehicle  buffet 
studies  (References  8  and  9)  has  been  applied  to  aircraft,  it  has  been  found 
that  the  fluctuating  buffet  pressures  on  a  highly  swept  airplane  wing  present 
very  different  characteristics.  In  general,  the  flow  pattern  on  the  upper 
surface  of  a  wing  dee^  in  buffet  does  not  exhibit  a  recognizable  flow  pattern, 
and  it  is  not  at  all  safe  to  assume  that  convected  flow  occurs  in  the  stream- 
wise  direction.  (See  Section  H  of  this  report). 

1.2  Buffet  Prediction  Method 

This  report  describes  the  development  of  a  method  to  predict  both  the 
fluctuating  pressures  and  the  structural  response  of  a  modern  fighter  air¬ 
plane  in  transonic  flight  when  it  penetrates  and  flies  beyond  the  buffet 
onset  boundary . 


1.2.1  Technical  Approach  -  The  variety  of  approaches  to  the  problem  of 
defining  the  input  loads  (i.e.,  the  spectral  and  spatial  energy  distribution 
of  the  fluctuating  pressures )ranges  from  purely  empirical  to  purely  theoretical. 
A  purely  empirical  approach  could  be  conceived  as  a  test  of  many  wing  geometries 
and  planfoi-ms  to  obtain  a  catalog  of  pressure  distributions  that  could  be 
applied  to  different  aircraft  wings.  The  main  disadvantage  of  this  approach 
is  the  inordinately  large  amount  of  wind  tunnel  testing  and  data  reduction 
that  would  be  required  to  provide  for  the  variety  of  wings  currently  in 
use  on  high  performance  fighters.  At  the  other  extreme,  a  purely  theoretical 
treatment  is  not  possib  le  because  the  state-of-the-art  in  theoretical 
aerodynamics  and  fluiu  mechanics  is  not  sufficiently  advanced  to  predict 
the  regions  and  levels  of  fluctuating  pressure  on  a  three  dimensional  wing. 

In  both  the  above  cases,  verification  of  predicted  buffet  response  loads 
would  remain  a  problem.  The  upproueh  selected  for  this  effort  uses  a 
combination  of  experimental  and  theoretical  inputs. 

Tliu  calculated  response, whether  it  be  wing  tip  acceleration  or 
.•root  sending  moment,  requires  an  input  excitation  force  (forcing  function) 
fud  a  structural  response  function.  The  development  of  the  prediction 
method  described  herein  is  bases  on  the  Interaction  effects  shown  in  Figure  1. 
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Ill  this  method  the  steps  are  defined  as  follows: 

(1)  First,  establish  the  nature  of  the  unsteady  pressures  acting  on 
the  wing. 

(2)  Second,  determine  the  surfaces  to  which  the  pressures  of  (l)  are 
to  be  applied. 

(3)  Next,  define  the  dynamic  characteristics  of  the  airplane  system. 

(M  The  fourth  step  provides  for  the  combination  of  pressure  -  area  - 

displacement  to  yield  a  forcing  function  in  the  form  of  general¬ 
ized  (modal)  buffeting  forces. 

(5)  The  forcing  function  is  then  applied  to  the  elastic  system  and 
the  required  buffet  response  loads  and  displacements  are  computed 
Step  (l),  of  course,  is  the  key  analytic  ingredient  which  is  obtained  by  a 
theoretical/experimental  method.  Step  (2)  is  based  on  wind  tunnel  flow 
visualization  studies,  and  step  (3)  may  be  done  either  analytically  or 
through  ground  vibration  tests.  Steps  (h)  and  (5)  are  computational 
procedures  that  have  been  obtained  by  modifying  an  existing  program.  The 
overall  objective  of  the  plan  is  the  development  of  n  first  generation 
buffet  prediction  method  for  use  in  preliminary  design.  The  method  offers 
capability  for  predicting  the  nature  of  the  unsteady  pressures  that  cause 
the  input  buffet  loads  as  well  as  the  response  leading  in  the  aircraft 
structure.  The  method  used  to  calculate  structural  response  loading  from 
input  buffet  dynamic  pressures  is  a  slight  modification  of  the  previously 
developed  method  of  Reference  10. 


SECTION  2 

:  TEST  PROGRAM 

The  experimental  phase  of  the  program  was  structured  around  two  series 
of  wind  tunnel  tests  that  were  conducted  on  a  10%  scale  YF-AE  airplane  model 
in  the  l6T  tunnel  at  AEDC,  The  purpose  of  the  first  test  (Series  I)  Was 
to  obtain  flow  visualization  data  on  the  upper  surface  of  the  wings  to  aid 
in  understanding  the  flow  patterns  and  to  aid  in  finalizing  locations  for 
fluctuating  pressure  transducers  (microphones).  The  purpose  of  the  second 
test  (Series  II)  was  to  obtain  measurements  of  the  time  histories  of  pressures, 
accelerations,  strains  and  forces  obtained  during  tunnel  riuns . 

2.1  Model  Design 

The  model  configuration  was  determined  by  the  YF-tE  version  of  the  basic 
F-tE  airplane  that  was  used  in  the  flight  test  program  described  in 
Reference  2.  A  full  span  10%  model  was  selected  in  order  to  make  optimum 
use  of  the  AEDC  l6l  tunnel  without  introducing  either  blockage  or  reflected 
shock  effects.  A  Task  Mk  XXI-A  balance  was  installed  between  the  model  and 
the  sting  to  measure  static  and  dynamic  components  of  forces  and  moments 
acting  on  the  model.  Figure  2  presents  a  three-view  of  the  model  showing 
overall  dimensions.  Both  the  inner  and  outer  wing  sections  (separated  by  a 
chordwise  line  at  the  snag  span  location)  are  equipped  with  leading  edge 
flaps.  Trailing  edge  flaps  are  located,  on  the  inner  wing  adjacent  to  the 
fuselage.  Rather  than  having  adjustable  flaps,  however,  the  model  was  fitted 
with  removable  flaps  of  different  angular  deflection  -  two  sets  for  the  leading 
edge  (0°  and  8°)  and  three  sets  for  the  trailing  edge  (0°,  T.5°,  and  15>°). 
b  •:  model  configurations  were  defined  on  the  . basis  of  leading  and  trailing 
edge  flap  combinations  as  shown  in  Table  1.  In  addition,  the  left-hand 
wing  was  fitted  with  a  row  of  oil  extrusion  holes  1/2  inch  apart  located 
along  a  line  that  went  from  about  1?  percent  chord  at  the  wing  root  to 
about  30  percent  chord  t  the  wing  tip.  These  holes  were  fed  from  a 
common  manifold  Imbedded  in  the  wing. 

2.2  Flow  Visualization  Tents 

The  first  phase  of  wind  tunnel  testing,  designated  8eri.es  I,  consisted 
>■-  -  obtaining  still  and  moving  pictures  of  oi  1  .flow  and- tuft  patterns  to 
identity  flow  characteristics  on  the  wing  and  to  aid  in  locating  fluctuating 
pressure  transducers.  The  left-hand  wing  was  equipped  with  the  oil 
orifices  and  had  its  upper  surface  painted  white  to  provide  contrast  with.. 


the  colored  oil.  The  upper  surface  of  the  right-hand  wing  was  painted  black 
to  provide  contrast  with  the  white  nylon  tufts  that  were  attached  in  a 
rectangular  grid  pattern  in  rows  approximately  one  inch  apart. 

The  flow  visualization  program  was  to  he  conducted  at  a  Reynolds  number 
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Re  =  3.75  x  10  per  foot  over  a  Mach  range  of  0.7  to  1.2  at  angles  of 
attack  between  -4°  and  20°.  However,  at  this  Reynolds  number  excessive 
dynamic  loads  were  experienced  at  ct^  greater  than  l4°tol6°.  To  permit 
reaching  the  higher.'  angles  of  attack  much  of  the  testing  was  conducted  at 
a  reduced  Reynolds  number  Re  =  1.25  x  10  ^  per  foot.  (A  limited  number  of 
runs  at  Re  =  2.5  x  10^  per  foot  were  also  made.)  Table  II  shows  the  test 
conditions  for  which  flow  data  were  obtained  (except  that  the  higher  angles 
of  attack  were  not  achieved  simultaneously  with  the  higher  Reynolds  numbers) 
A  photograph  showing  both  oil  flow  and  tufts  during  tunnel  flow  is  shown 
in  Figure  3.  Additional  test  details,,  including  aerodynamic  force  data, 
are  given  in  References  11  and  12. 

2.3  Ilynamic  Instrumentation  on  Model 

After  the  flow  visualization  test  and  prior  to  model  vibration  tests, 
the  dynamic  instruments  were  installed  on  the  model.  This  instrumentation 
consisted  of  fluctuating  pressure  transducers  (microphones),  accelerometers 
and  strain  gages.  The  microphones  were  installed  on  the  upper  surface  of 
the -wing  and  tail  to  measure  fluctuating  surface  pressures.  Three  strain 
gages  were  installed  on  the  wing  and  one  was  installed  on  the  horizontal 
tail  to  measure  bending  and  torsional  strains  induced  during  buffet.  In 
addition  to  an  accelerometer  located  at  the  model  c.g, ,  a  pair  of  high 
frequency  accelerometers  was  located  near  the  wing  tip  to  measure  bending 
and  torsion  acceleration  response  to  buffet  pressure.  Figure  4  shows  the 
overall  model  instrumentation  while  Figure  5  shows  the  detailed  location 
of  the  sensors  on  the  wing.  Ail  wing  mounted  transducers  were  recessed  in 
the  wing  to  preserve  a  smooth  wing  contour.  Leads  were  routed  through 
grooves  In  the  wing,  into  the  fuselage  and  hence  into  the  sting. 

A  description  of  the  date  acquisition  system  used  by-  AE1KJ  at  the  loT 
tunnel  is  given  in  paragraph  3.1  of  the  next  section. 

2.3,1  Denser  Tioncvi pt iun  -  The  dynamic  sensors  used  on  the  101#  VF-4E 
model  are  described  us  follows; 

Mi  crophones  - 

Disc  mounted  -  Kulite  \'/U  CQL-060-25 


Flush  mounted  -  Kulite  P/N  Uiii-150-25 
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Table  I  -  Model  Configurations 


Configuration 

Leading  Edge  Flap  Deflection 
(degrees) 

1  (Clean) 

0  - 

2 

8 

3 

8 

Trailing  Edge  Flap  Deflection 
(degrees) 


5 

0 

16 

6 

8 

15 

Table 31  -  Test  Conditions  for  Flow  Visualization  (Series  I) 


Configuration  M 


0,7, 0.9, 1.0, 1,2 


0.7, 0,9, 1,0, 1.1, 1,2 


07,0.9,1.0.1.2 


07,0.9, 1.0, 1.2 


07,0.9. 1.0, 1.2 


07. 0,85. 0.9, 1.0. 1.2 


5tw  l«»ule  1,  Uolururuiu  11  fur  Additional  DetuHs 


Range of  ay 
(deg) 


-  2,  -»  20 


-2,  +?0 


..~2,  +20 


-4,  +20 


-3,  +20 


Re/ft  (xlQ*6) 


1.25,375 


1.20. 2.5, 375 


1.26,375 


1.25,375 


1.26, 376 


1.25, 375 


k 


«-y 


Accelerometers  - 
Endevco  2222k 
Strain  Gages  - 


Two-element  rosettes  (2  per  Wheatstone  bridge  to  measure  torsion) 
Bean  BAB-06-125XS-350 

Uniaxial  gages  (1+  per  Wheatstone  bridge  to  measure  bending) 
Micro-Measurements  EA-06-090DH-350 
Micro-Measurements  EA-06-090EG-350 

2‘3*2  Microphone  Grid  -  A  grid  of  twenty-eight  wide  band  fluctuating 
pressure  transducers  was  installed  on  the  model  prior  to  the  buffet  data 
acquisition  tests.  The  grid  was  based  on  the  following  criteria: 

(1)  Microphones  would  be  densely  grouped  in  areas  of  high  turbulence 
and  sparsely  grouped  in  regions  of  low  turbulence. 

(2)  The  microphones  would  be  arranged  so  that  key  flow  characteristics 
.could  be  readily  determined. 

(3)  Some  microphones  would  be  located  to  provide  direct  comparison 
with  flight  test  data. 

(M  Fluctuating  pressure  on  the  tail  should  be  monitored. 

The  microphone  locations  shown  in  Figures  4  and  5  were  based  on  these 
criteria  and  the  results  of  flow  visualization  tests.  The  forward  most 
locations  on  the  inner  wing  panel  could  not  extend  forward  of  the  12% 
chord  line  because  of  interference  with  the  flap,  and  those  on  the  outer 
wing  panel  could  not  be  forward  of  about  the  30$  chord  line  because  of  thin¬ 
ness  of  the  wing  and  interference  with  the  oil  flow  manifold.  Also,  the 
locations  of  transducers  17  and  10  could  be  located  no  further  aft  than 
the  90$  chord  (because  of  wing  thinness)  although  it  was  desirable  to  corre¬ 
late  the  outputs  with  flight  test  pickups  located  at  95$  chord.  However, 
transducers  17  and  18  do  coincide  with  those  of  flight  test  pickups  on  the 
90$  chord.  The  functional  groupings  of  the  sensors  are  listed  in  Table  III 
together  with  the  purpose  each  grouping  war  intended  to  serve. 

The  transducers  in  Group  I  formed  four  ehordwise  rows.  Three  of  these 
rows  were  about  equally  spaced  on  the  outboard  half  of  the  wing.  Emphasis 
•  /'Ivan  to  the  outboard  half  because  os  angle  of  attack  is  increased  this 
region  is  the  first  to  have  separated  airflow.  The  location  of  the  inboard 
row  consisting  of  5,  14,  and  21  was  fixed  by  pressure  transducers  on  the 
flight  test  airplane. 


Table  HI 

Functional  Grouping  for  Fluctuating  Pressure  Transducers 


Group 

Transducers 
(See  Figure  4) 

Purpose 

1 

I, 2,4,6,6,7,9,10, 

II, 13,14,16,18, 
20,21 

Coarse  Grid  for  Defining  Major  Flow 
Features  on  Wing 

2 

11,15,16,17,18 

21,26,27,28 

21,22,23,24,25 

Fine  Grid  for  Defining  Flow  Features 
on  Outer  Wing 

Wing-Tail  Correlation 

Close  Correlation  on  Flap  and  Check 
of  Pressure  Distribution 

3 

2,3,4,7,8,9,11,12, 

13,18,19,20 

Investigate  Localized  Effect  of  Snag 

4 

17,18,21,24 

Compare  with  Flight  Test  Data 

The  first  grid  in  Group  2  was  to  examine  the  streamwise  and  spanwise 
correlation  between  closely  spaced  points  on  the  outer  wing  panel.  The 
second  grid  in  Group  2  was  for  determining  correlations  between  fluctuating 
pressures.-on  the  wing  and  tail,  and  between  those  on  the  two  sides  of  the 
tail.  The  third  grid  was  to  provide  a  check  of  the  mathematical  model  for 
the  fluctuating  pressures  and  contained  five  transducers  on  the  trailing 
edge  flaps.  This  check  could  not  be  made  because  a  strain  gage  on  the  flap 
support,  needed  for  the  correlation,  malfunctioned. 

The  purpose  of  the  transducers  in  Group  3  was  to  investigate  the 
localized  effect  of  the  snag. 

The  locations  of  the  fluctuating  pressure  transducers  in  Group  k  were 
the  same  as  those  in  the  flight  test  airplane  (Reference  2). 

Twenty- four  of  the  twenty-eight  microphones  used  on  the  model  were 
installed  in  mounting  discs  (see  Figure  lH  of  Reference  13)  that  recessed 
the  diaphragms  and  protected  them  from  any  possible  direct  stream  impinge¬ 
ment.  Comparative  frequency  calibrations  showed  that  frequency  response 
was  not  affected  by  the  mounting  discs.  The  calibration  is  discussed  in 
Reference  13.  No  microphones  were  lost  due  to  diaphragm  damage  during  the 
Series  II  tests. 

2.H  Vibration  Tests 

A  vibration  test  was  made  after  installation  of  sensors,  to  determine 
the  vibration  modes  and  frequencies  that  might  influence  buffet  data  measured 

on  the  model.  It  was  believed  that  wide-band  buffet  pressure  excitation 
could  excite  and  amplify  resonant  structural  modes  that  would  in  turn  in¬ 
duce  pressure  peaks.  Such  peaks  should  obviously  not  be  included  in  mathe¬ 
matical  modeling  of  the  pressure.  The  test  was  done  in  two  parts  at  two 
different  locations  -  the  first  part  was  conducted  in  the  MCAIR  laboratories 
and  consisted  of  a  comprehensive  evaluation  of  model  modeshapes  and  frequencies. 
These  data  were  obtained  with  the  setup  shown  in  Figure  6.  The  model  was 
mounted  through  the  model  balance  flexures  to  a  rigid  sting,  and  excitation 
was  provided  by  a  pair  of  electromagnetic  shakers.  By  properly  phasing  the 
shaker  force  outputs  (i.e,,  in-phase  or  lOO°  out-of-phase ) ,  either  symme¬ 
trical  or  anti -symmetrical  vibration  modes  could  be  excited.  The  teat 
procedure  consisted  of  (l)  sweeping  through  a  range  of  frequencies  between  0 
and  1500  Hz  and  recording  the  output  of  monitor  accelerometers,  (2)  identify¬ 
ing  the  modes  represented  by  the  peaks  of  the  swept  response  and  (3)  mapping 
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dominant  modes.  Table  IV  lists  and  identifies  model  resonant  frequencies. 

The  remaining  part  of  the  vibration  test  was  made  at  AEDC  with  the 
model  mounted  through  the  balance  flexures  to  the  wind  tunnel  sting.  The 
vibration  test  setup  shown  in  Figure  7  was  used  to  determine  the  sting/model 
vibration  characteristics.  Sting/model  frequencies  measured  in  this  test 
are  listed  in  Table  V.  Details  of  the  vibration  tests  are  presented  in 
Reference  14. 

2.5  Buffet  Data  Acquisition  Tests 

The  Series  II  buffet  data  acquisition  tests  consisted  of  a  series  of 
tunnel  runs  at  various  Mach  numbers  and  angles  of  attack  to  gather  pressure , 
acceleration,  strain }and  force  and  moment  coefficient  data  related  to  air¬ 
craft  buffet. 

The  test  was  conducted  at  Mach  numbers  between  0.7  and  1.10  and  most  of 
the  data  were  obtained  at  a  Reynolds  number  of  3.75  millions  per  foot. 

The  angle-of-attack  range  was  from  -1.5  to  16  degrees  at  zero  degrees  yaw. 

Data  also  were  taken  with  configuration  1  (clean)  through  a  Reynolds 
number  range  of  1.25  through  5.25  millions  per  foot  at  several  angles  of 
attack.  A  summary  of  parameter  variations  is  presented  in  Table  VI. 

The  general  operating  procedure  was  to  record  simultaneously  60  seconds 
of  VIDAR  analog  data,  and  then  record  tunnel  and  model  static  pressure  data 
digitally.  This  sequence  was  done  at  an  angle  corresponding  to  zero  lift 
(Cr  =  0)  and  at  a  number  of  angles  in  the  region  of  buffet  onset  and  beyond. 
At  some  additional  angles  of  attack,  digital  data  only  were  recorded  for  the 
purpose  of  better  definition  of  overall  aircraft  aerodynamic  curves,  A 
complete  list  of  the  angles  at  which  data  were  recorded  is  presented  in 
Table  2  of  Reference  13. 
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SECTION  3  .  -.v;V:  7V 

DATA  ACQUISITION  AND  REDUCTION  ' 

This  section  describes  the  data  acquisition  and  data  reduction  tech¬ 
niques  utilized  in  obtaining  the  power  spectral  density  plots  that  are 
employed  in  deriving  a  mathematical  model  of  the  fluctuating  pressure.  The 

i 

raw  data  were  obtained  during  the  Series  II  wind  tunnel  tests  conducted  at 
the  AEDC  in  the  l6T  transonic  tunnel. 

3.1  On-Line  Data  Acquisition  and  Reduction  System  at  AEDC 

The  l6T  tunnel  was  equipped  with  a  readout  system  for  obtaining 
and  digitizing  the  mean  and  rms  outputs  of  the  strain-gage  and  piezoelectric 
type  transducers  installed  in  the  model/sting  balance  assembly.  On-line 
data  consisted  of  a  tabulation  of  the  mean  (static)  and  rms  (dynamic)  response 
values  in  engineering  units  of  the  transducer  signals “during  each  tunnel  run 
performed  at  a  specific  Mach  number,  model  angle  of  attack,  and  model  config¬ 
uration. 


The  on-line  data  acquisition  system  provided  FM  tape  records  of  the  raw 
transducer  outputs  for  each  wind  tunnel  run  for  off-line  data  analysis.  The 
analog  data  were  recorded  on  a  multiplex  tape  recording  system  utilizing  a  CEC 
Model  VR  TOO  wide-band  tape  recorder  transport,  Vidar  multiplexing  electronics, 
and  ill  track  FM  magnetic  tapes.  A  tape  recording  speed  of  60  inches  per  second 
(ips)  with  discriminator  center  frequencies  of  l6o,  220,  280,  340,  400,  and 
460  KHz  provided  effective  bandwidths  of  8  KHz.  A  description  of  the  on-line 
data  acquisition  and  reduction  system  used  for  the  l6rf  tunnel  is  presented  in 
Reference  15* 

3 . 2  Data  Reduction  at  MCAIR 

Off-line  analysis  of  the  FM  magnetic  tape  records  of  the  dynamic 
data  recorded  during  the  buffet  dynamic  loads  program  was  performed  at  MCAIR, 
Since  the  raw  data  tapes  were  recorded  by  a  multiplex  frequency , modulated 
recorder,  the  tapes  were  played  back  through  discriminators  to  separate  the 
individual  data  signals.  The  individual  signals  were  then  stripped  out  on 
a  recording  oscillograph  and  time  histories  were  obtained  for  all  outputs  of 
the  acceleration,  strain,  and  dynamic  pressure  transducers.  The  oscillograph 
.  ,c.(  ;  were  subsequently  edited  for  data  validity  and  thirty-five  transducer 
pairs  were  chosen  to  be  analyzed.  Choice  of  the  transducer  pairs  was  guided 
by  the  kinds  of  information  being  sought:  the  spectral  content,  convection,  • 
and  spatial  decay  of  fluctuating  pressures  acting  on  the  wing  and  the  spectral 


Ov>ntent  of  and  correlation  between  various  model  responses .  The  data  were 
-reduced  for  the  1 6  test  conditions  that  are  outlined  in  Table  VII.  Most  of 
these  test  conditions  were  chosen  to  exhibit  the  effect  of  Mach  number,  angle 
of  attack,  and  model  configuration  on  buffet.  Other  conditions  were  chosen 
to  match  those  of  flight  buffet  tests,. or  to  show  the, nature  of  background 
noise  (at  CT  =  0), 

3.2.1  Oscillograph  Traces  -  Time  histories  of  typical  dynamic  pressures, 
accelerations,  and  bending  moments  recorded  during  the  test  are  presented  in 
Figures  8  through  11.  Since  the  fitting  of  these- traces  by  deterministic 
functions  seems  unlikely,  the  data  were  treated  as  random.  Moreover,  the  data 
were  taken  to  be  stationary  since  no  changes  in  the  statistical  characteris¬ 
tics.-.  are  obvious  from  visual  examination  of  the  traces. 

i  3.2.2  Broad  Band  FM  Recording  -  After  editing,  the  data  selected  for 
reduction  were  re-recorded  on  broad  band  FM  tape.  Since  the  center  frequencies 
and  effective  bandwidths  of  the  MCAIR  discriminators  were  twice  those  used  to 
record  the  data  at  AEDC,  the  multiplexed  data  were  played  back  through  the 
-MCAIR  discriminators  at  120  inches  per  second.  The  signals  obtained  from  the 
discriminators  were  recorded  on  broad  band  FM  tape  at  30  ips  and  a  center 
frequency  of  27  KHz  in  a  5  KHz  frequency  band. 

The  actual  data  reduction  was  performed  on  tape  loops  obtained  from  the 
30  ips,  27  KHz  FM  tapes.  These  loops,  made  to  minimize  phase  shifts  and  to 
insure  coherence  of  the  transducers  in  a  given  pair,  were  used  for  all  spectral 
analysis  at  MCAIR.  They  were  played  back  at  15  ips  with  a  center  frequency  of 
13.5  KHz.  Thus,  the  effective  frequency  response  of  the  data  analysis  was 
from  0  to  2.5  KHz. 

3.3  Mathematical  Functions,  Spectral  Analysis 

As  discussed  in  Section  3.2.1,  the  dynamic  data  studied  in  this  program 
appear  to  be  random  and  stationary.  The  analysis  of  such  data  is  generally 
given  in  terms  of  correlation  and  spectral  density  functions,  which  are 
Fourier  transforms  of  one  another.  These  functions  are  discussed  below, 
and  further  details  may  be  found  in  References  16  and  17. 

3.3.1  Correlation  Func lions  -  The  correlation  functions  for  random 
data  are  defined  as  follows: 

Autocorrelation  function 
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Test  Conditions  for  Reduced  Dynamic  Data 


Model 

Configuration 

Mach 

Number 

Reynolds 
Number 
(per  foot) 

Wing 

Angle-of-Attack 

(deg) 

1 

0.7 

3.75  x  106 

-  0.85 

1 

0.7 

3.75  x  106 

12.08 

1 

0.85 

3.75  x  106 

10.02 

1 

0.9 

3.75  x  106 

7.98 

1 

0.9 

3.75  x  106 

12.08 

1 

0.9 

3.75  x  106 

14.08 

1 

0.9 

3.75  x  106 

-  0.85 

1 

0.8 

3.75  x106 

8.99 

1 

0.8 

3.75  x  106 

9.95 

1 

1.0 

3.75  x  106 

0 

1 

1.0 

3.75  x  106 

12.0 

5 

0.7 

3.75  x  106 

11.0 

5 

0.9 

3.75  x  106 

11.0 

5 

0.9 

3.75  x  106 

14.0 

G 

0.7 

3.75  X  106 

11.0 

6 

0.8 

3.75  x  IQ6 

9.0 

Cross-correlation  function  . 

12  -i  T 

R  (t)  =  Lim  ~~  /  f  (t)  f  (t+x)dt  (3-2) 

T  -*  co  "  x  d 

where 

f^(t)  =  time  function  1 

fg(t)  =  time  function  2 

T  =  one-half  of  integration  period 

R(x)  =  correlation  function  as  function  of  time  delay  x. 

Autocorrelation  and  cross-correlation  functions  were  computed  using  a 
Princeton  Applied  Research  analyzer.  The  analyzer  performed  the"  following 
function 

T 

R^U)  =  Lira  /  f  (t)f?(t  -  x)dt  (3-3) 

T  -►  oo  i  d 


The  primary  difference  between  this  function  and  the  one  expressed  in 
Equation  (3-2)  is  the  time  delay  -x.  A  negative  time  delay  does  not  affect 
the  autocorrelation.  The  cross-correlation  is  affected  as  follows 

bis(0  -a12’  (-,) 

(34) 


The  cross-correlation  function  can  he  changed  to  correspond  to  the  format 

of  Equation  (3-2)  by  changing  the  sign  of  the  abscissa. 

3.3.2  Power  and  Cross  Power  Spectral  Pen,  tty  Functions  -  The  power 

spectrum  of  a  stationary  random  function  is  the  Fourier  transform  of  the 
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correlation  function. 

Power  Spectral  Density 
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Sn(f>  =  /  RU(T)o-ia”fT  <1,  (3-5) 

Cross  Power  Special  Density 
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where 

S11(f)  =  power  density  spectrum  as  a  function  of  cyclic  frequency  (f) 
defined  for  -»<  f<»  and  with  units  (amplitude  l)2/Hz 
S1^(f)  =  cross  power  density  spectrum  as  a  function  of  (f)  defined 
for  -«<  f<»  and  with  units  (amplitude  1) (amplitude  2)/Hz 

The  corresponding  inverse  transforms  are 


Rax(t)  =  /  Si:L(f)ei2 'Tf  df 


(3-7) 


B12(t)  »  /  S12(f)ei2lTTf  df 


(3-8) 


The  mean  squared  value  of  a  time  function  is  the  value  of  the 
correlation  function  at  t  =  0,  This  function  can  be  expressed  as  the 
integrated  power  spectrum 

1  CO 

of,  =  Fn(0)  =  Lim  —  !  f?(t)  dt  a  /  SU(f)  df 

11  t  +  »  -T  1  _oq 

0Q  W 

=  2  /  SU(f)  df  a  /  Gn(f)  df  (3-9) 


where 

G^( f)  =  the  "physical”  one-sided  power  spectral  density  defined  for  0  <f. 
The  above  equation  is  used  in  the  physical  determination  of  the  power 
spectrum  and  the  cress  power  spectrum. 

The  analog  technique  used  to  compute  power  spectral  density  is  based 
on  the  following  operations:  filtering,  true  square  detection, and  integration. 
The  combination  of  these  three  functions  allows  the  computation  of  power 
spectral  density  for  random  funetipns.  A  block  diagram  of  the  operation  of 
the  analog  analyzer  follows : 
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If  the  time  function  is  filtered  and  the  bandwidth  of  the  filter  is 
small,  the  power  spectral  density  within  the  bandwidth  can  be  assumed  to 
be  constant.  Equation  (3-9)  can  then  be  rewritten  as  follows: 

Af 

c  2  =R11(0)=  Lim  S  fF2(t)dt  =  G^(f)  /  df  =  Gn(f)Af  (3-10) 

11  t  ->  “  -T  ~Af 

2 

where 

G^(f)  =  power  spectral  density  of  time  function  1  at  frequency  f 
Af  =  filter  bandwidth,  Hz 

The  power  spectral  density  can  now  be  expressed  in  terms  of  the 


filtered  time  function 

Gn(f)  =  Lim 

Af**0,  T-*» 


2TAf  ^ 


(3-11) 


A  similar  relation  exists  for  the  cross  power  spectral  density.  The 
cross  power  spectrum  is  expressed  in  tenns  of  the  in-phase  (real.  Re) 
and  whe  quadrature  phase  (imaginary,  Im)  functions.  (These  functions  are 
also  called  the  co-spectrum  C  and  the  quad-spectrum  Q  respectively.) 


Re  [G12(f)l  «  Um  /  fxp(t)  f^t)  dt  «  C12(f) 

Af—Q,  T*®*  -T  / 


(3-12) 


imaginary 


Im  (G12(f)J  »  Lim  /  f\»(t)-fw(t)dt  «  Q12(t) 

Afu-O/f**  ■■  -T  ■  **  , 


(3-13) 


where 

flp(t)  »  filtered  time  function  f,(t) 
i%y(t)  «  filtered  time  function  iV,(t) 

.  «  filtered  time  function  of  f^U)  phase  shifted  90° 
f,i%)  «*  the  ‘'phyulcui”  one-sided  cross  power  spectral  density  defined 
.  for  0‘  f . 

IM 

*ihe  ar^ia  power  spectral  density  can  then  he  expressed  as  the  following 
complex  quantity 

~  *  -  1  la 

*  3 1?t f 5  -  i  Qw(t)  ,  (3-lfc) 


or,  witten  in  terms  of  magnitude  and  phase, 
G12(f)  =  I  Gi2(f)  I  e  12  '•  •  • 


(3-15) 


where 


l^f)  |=[[Cl2(f)-]2 .+  [Q12(f)]2]l/2 


012^f)  =  tan-1  [Q  (f)/C  (f)] 

and  i  =  J-l 


(3-16) 


(3-17) 


Figure  12  shows  the  data  flow  and  steps  involved  in  reducing  raw  data 
to  obi  in  PSD  and  CPSD  plots.  Power  spectral  density  (PSD)  analysis  was 
performed  on  a  Honeywell  Model  ASA9300B  Automatic  Spectrum  Analyzer  which 
can  analyze  a  pair  of  transducers  simultaneously.  This  analyzer  consisted 
of  a  programmable  tuning  oscillator,  a  pair  of  filters  having  four  possible 
bandwidths  (2  Hz,  5  Hz,  10  Hz  and  50  Hz),  a  true  square  detector,  an 
integrator,  and  a  logarithmic  converter.  The  tuning  oscillator  of  the 
Honeywell  system  was  programmed  to  step  through  a  spectrum  from  0  to  2000  Hz 
utilizing  bandwidths  and  integration  times  as  follows: 


Frequency 
0  Hz  -  30  Hz 
30  Hz  -  150  Hz 
150  Hz  -  800  Hz 
800  Hz  -  2000  HZ 


Filter  Bandwidth 


10  Hz 


50  Hz 


Integration  Time 
20.0  sec. 
20.0  sec. 
10.0  sec. 

3.0  see. 


The  filtered  data  were  then  further  processed  to  obtain  the  PSD.  The 
programmable  oscillator  also  provided  a  D.C,  voltage  signal  proportional 
to  the  center  frequency  of  the  bandwidth  being  analyzed.  The  PSD’s,  along 
with  the  frequency  related  voltage,  were  recorded  digitally  on  a  DYMEC 
Model  DY-201QH  Data  Acquisition  System,  which  was  compatible  with  an  IBM 
360/85  input  format. 

A  parallel  .tap  was  made  between  the  narrow  bandpass  filter  and  the 
true  square  detector  in  the  Honeywell  system  to  obtain  filtered  data  in 
each  stepped  bandwidth  for  coincident  and  quadrature  analysis.  The  coin¬ 
cident  and  quadrature  analysis  was  performed  using  Spectral  Dynamics  Model 
3D  1012  Tracking  Filter,  a  Model  SD  1010  Carrier  Generator,  and  a  Model  3D 
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109A  Co/Quad  Analyzer.  The  center  frequency  of  the  carrier  generator  was 
referenced  to  a  7.5  KHz  signal  used  in  the  Honeywell  system.  The  tracking 
filter  (equipped  with  100  Hz  bandwidth  filters)  was  used  to,  provide  com¬ 
patibility  between  the  spectrum  analyzer,  operating  with  a  carrier  frequency 
of  7.5  KHz,  and  the  co/quad  analyzer  having  a  carrier  frequency  of  100  KHz. 
The  100  Hz  filters  caused  no  additional  filtering  of  the  data  since  they  were 
twice  the  possible  programmed  bandwidth  of  the  Honeywell  system.  Thus,  a 
coincident  and  a  quadrature  value  were  obtained  for  each  bandwidth  ’’stepped 
off"  by  the  Honeywell  system.  The  values  of  the  coincidence  and  the  quad¬ 
rature  were  recorded  digitally  utilizing  the  DYMEC  Data  Acquisition  System. 
The  DYMEC  was  also  able  to  record  a  ma..  :1  input  channel  which  provided 
the  recording  of  run  number,  transducer  pairs,  system  gain  factors,  band¬ 
width,  and  integration  times. 

The  IBM  360/85  computer  was  used  to  read  the  digital  records,  convert 
all  data  to  engineering  units,  and  to  <'mpute  the  magnitude  and  the  phase  of 
the  cross  power  .spectial  density,  the  square  root  of  thf  coherence  function, 
and  the  cumulative  root  mean  square  values  through  the  PSD  spectrum  for 
both  transducers.  The  computer  presented  the  data  in  two  formats.  The  first 
was  in  the  form  of  a  computer  printout  for  each  "stepped  off"  bandwidth 
throughout  the  0  Hz  to  ?00Q  Hz  frequency  spectrum  of  each  parameter  calcu¬ 
lated.  The  second  format  consisted  of  a  digital  "history  tape"  for  later  use 
on  the  Model  663  Cal  Comp  Plotter.  The  history  tape  provided  the  means  by 
which  data  could  be  saved  for  later  recall  and  plotting,  if  necessary. 

3 . 4  Accuracy  of  Data 

Uncertainties  in  the  dynamic  data  analyzed  in  this  test  program  result 
from  (l)  the  instrumentation  errors  associated  with  the  hardware  that  was 
uaod  to  convert  recorded  raw  data  into  power  spectra  and  (2)  statistical 
inaccuracy.  instrumentation  accuracy  is  estimated  to  be  +8$  of  the  data 
presented.  The  larger  uncertainly  is  due  to  nfcatiatieal  inaccuracy,  which 
is  presented  for  the  PSD  and  CPSD  data  in  Tables  VilX(a)  and  Vin(b) 
respectively.  As  indicated  in  this  table,  the  standard  deviation  (error) 
introduces  uncertainty  in  the  cross  power  spectral  data  as  large  as  22.1*$. 


Table  Sin 

Normalized  Standard  Errors 


(a)  Power  Spectral  Density 


Spectrum 
Section  (Hz) 

Bandwidth 

(Hz) 

integration 
Time  (sec) 

Statistical 
Degrees  of 
Freedom* 

Normalized 

Standard 

Error** 

1-30 

2 

20 

80 

15.8% 

30-150 

5 

20 

200 

10.0% 

150  -  800 

10 

10 

200 

10.0% 

800  -  2000 

50 

3 

300 

8.2% 

‘Degrees  of  Freedom  =  2  (Bandwidth)  (Integration  Time) 

‘‘Normalized  Standard  Error,  %  =  100/7  (Bandwidth)  (Integration  Time) 
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(b)  Cross  Power  Spectral  Density 


Spectrum 
Section  (Hz) 

Bandwidth 

(Hz) 

Integration 
Time  (sec) 

Statistical 
Degrees  of 
Freedom 

Normalized 

Standard 

Error 

1  -  30 

2 

10 

40 

22.4% 

30-160 

6 

10 

100 

14.1% 

160-800 

10 

5 

100 

14.1% 

800  2000 

HO 

l.h 

150 

11.5% 

& 
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SECTION  V 
DATA  ANALYSIS 


This  section  presents  discussion  and  some  results  from  analysis  of  data 
obtained  during  the  flow  visualization  (Series  I)  and  buffet data  acquisition 
(Series  II )  wind  tunnel  tests.  The  first  part  of  the  analysis ,  related  to 
aerodynamic  flow,  had  as  its  principal  objectives:  (l)  determine  the  nature 
of  the  flow  on  the  upper  wing  as  it  varies  with  Mach  number,  with  angle  of  ' 
attack,  and  with  Reynolds  number,  (2)  assess  the  effect  on  pressure  distri-' 
bution  of  vortices,  (3)  define  buffet  onset  and  (4)  locate  fluctuating  pres¬ 
sure  transducers.  The  second  part  (data  from  Series  II  test),  dealing  with 
rms  and  PSD  data,  had  the  following  main  objectives:  (l)  define  the  trends 
ACprms  (overall  rms  value  of  fluctuating  pressure)  with  angle  of  attack 
and  Mach  number,  and  with  position  on  the  wing,  (2)  determine  the  distribu¬ 
tion  of  fluctuating  pressures  with  respect  to  frequency  and  position  < n  the 
wing  during  changes  in  angle  of  attack  and  Mach  number,  and  (3)  make  a  com¬ 
parison  of  wind  tunnel  and  flight  data. 

4.1  Flow  Analysis 

The  flow  characteristics  on  the  upper  wing  surfaces  were  monitored 
visually  by  tufts  and  oil  flow  to  determine  flow  patterns  and  aid  in  locating 
transducers.  Data  reduction  was  done  using  color  movies  and  black  and  white 
still  photographs.  A  summary  of  the  type  of  information  obtained  is  presented 
below  and  details  of  flow  analysis  are  given  in  Appendix  A.  Additional  still 
photographs  are  given  in  Reference  12. 

4.1.1  Tufts  -  Nylon  tufts  were  attached  to  the  upper  surface  of  the 
right  wing  in  rectangular  rows  1  inch  apart.  For  analysis,  tuft  motion  was 
categorized  according  to  degree  of  turbulence  based  upon  tuft  motion  and 
direction. 

4.1.2  Oil  Flow  -  Oil.  flow  photographs  were  analyzed  in  a  similar 
manner  although  the  eritieria  for  Judging  degree  of  turbulence  (and  hence 
separation)  were  different .  Flow  where  oil  streaks  were  largely  oriented 
sitreamwiae  and  where  most  of  the  oil  was  rapidly  swept  away  was  Judged  to 
ixs  attached.  At  the  opposite  extreme,  highly  turbulent  (separated)  flow  is 
Judged  to  occur  where  oil  collects  and  appears  to  have  a  poorly  defined  flow 
direction.  Separation  boundaries  were  identified  by  using  both  tufts  and 


oil  flow. 
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4.1.3  Flow  Characteristics  -  The  flow  on  the  F-4  model  wing  is 
strongly  influenced  by  the  leading  edge  discontinuity  called  the  snag. 

Figure  13  is  an  oil  flow  photograph  made  at  M  =  0.9  and  aw  =  12°  that  shows 
both  the  tuft  and  oil  flow  patterns  for  a  representative  buffet  condition. 

In  general,  flow  separation  began  on  the  outer  wing  panel,  and  progressed 
inboard  with  increasing  angle  of  attack. 

4.2  Analysis  of  rms  Data 

4.2.1  Noise  in  Wind  Tunnel  -  It  is  known  that  a  wind  tunnel  produces 
noise.  Special  runs  were  made  for  the  clean  configuration  at  C.  =  0  to 
determine  the  nature  of  this  background  noise.  In  this  zero  lift  condition 
the  flow  is  attached  and  laminar  and  pressure  fluctuations  are  attributed 
to  perturbations  in  the  stream  flow.  Noise  level  measurements  are  compared 
with  turbulence  measurements  in  the  following  paragraphs ,  and  it  is  shown 
that  in  general  the  noise  level  is  not  high  compared  to  the  buffeting  pres¬ 
sure  level. 

4.2.2  Pressure  rms  Measurements  -  In  the  analysis  of  rms  pressure 

measurements  presented  in  this  report,  emphasis  was  placed  on  configuration 

1  for  examining  the  effect  of  Mach  number  and  angle  of  attack.  At  Mach  0.9 

data  were  analyzed  at  =  8° ,  12° ,  and  l4°  to  3tudy  the  effect  of  angle  of 

attack;  at  12°  wing  angle  of  attack,  data  were  analyzed  at  M  =  0.7,  0.9,  and 

1.0  to  study  the  effect  of  Mach  number.  For  configuration  5,  the  effects  of 

angle  of  attack  and  Mach  number  were  examined  by  analyzing  data  at  M  =  0.7, 

a  =  11°;  M  =  0.9,  a  =  11°:  and  M  =  0.9,  a  =  l4°.  To  provide  an  indication 
w  w  w 

of  the  effect  of  model  configuration,  data  were  analyzed  for  Configuration  6 
at  M  “  0.7,  a  =  11°  for  comparison  with  data  at  similar,  conditions  .for ' 
Configurations  1  and  5. 

Figure  l4  shows  the  variations  of  ACp  with  percent  chord  for  the 
model  at  CL  -  0  and  Mach  0.7.  This  value  (between  0.01  and  0.02)  is 
assumed  to  represent  tunnel  noise  and  will  serve  as  a  Mach  0, J  reference  for 
the  buffet  data.  Buffet  pressures  significantly  higher  than  this  will  be 
<>.; sen tULly  unaffected  by  tunnel  noise.  The  variation  of  ACprjng  at  Mach  0.7 

arid  a  »  32°  is  shown  in  Figure  35.  The  values  are  usually  well  above  the 

w  . 

level  and  below  0.1  except  for  the  aft  part  of  the  78  percent  semispan, 
v/lvre  they  are  above  0.2.  This  high  level  is  associated  with  the  snag  vortex 
an  discussed  in  Appendix  A. 
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Figure  1 6  shows  variation  of  ACp^mg  with  percent  span  at  M  =  0.9  and 
'  C^  =  0.  Except  for  one  point,  these  values,  are  below  0.02,  which  is  taken 
as  the  tunnel  noise  level  at  Mach  0.9.  Figures  17,  18,  and  19  are  plots  of 
ACp^ms  versus  percent  chord  at  Mach  0.9  and  wing  angles-of-attack  of  8°,  12° 
and  l4°,  respectively.  Most  of  the  ACp^ms  data  are  below  0.07.  The  excep¬ 
tions  are  the  78  percent  semispan  and  90  percent  chord  at  a  =8°  and  the 

w 

78  and  93  percent  semispan  at  a  =  12°  and  14°.  In  the  latter  . two  cases, 
ACprffis  gets  as  high  as  0.l6,  These  high  values  are  thought  to  be  associated 
with  the  snag  and  tip  vortices  as  discussed  in  Appendix  A. 

The  tunnel  noise  level  at  Mach  1  is  shown  by  the  plot  of  ACr,  versus 
percent  chord  in  Figure  20  to  be  about  0.01.  The  x'ew  values  above  this  level 
are  thought  to  be  due  to  snag  and  tip  vortices  and  some  flow  perturbation 
caused  by  the  engine  inlet.  Figure  21  shows  that  at  Mach  1  and  =  12°, 
ACprmg  i-s  below  0.0 6  at  all  the  measurement  locations.  It  thus  appears. that 
the  larger  buffeting  pressures  occur  below  Mach  1. 

Variations  of  ACp  g  over  the  wing  with  the  trailing  edge  flaps  deflected 
15  degrees  are  shown  in  Figures  22,  23,  and  24.  A  comparison  of  Figures  15 
and  22  shows  little  effect  due  to  trailing  edge  flap  deflection  at  Mach  0.7 
and  aw  in  the  neighborhood  of  12°.  Comparing  Figures  18  and  23  shows  that 
at  Mach  0.9  and  around  a  =  12° .deflecting  the  trailing  edge  flaps  15° 

•raises  ACp^^  roughly  30  percent  at  the  78  percent  semispan  and  inboard  and 
lowers  it  slightly  at  the  tip.  At  Mach  0.9  and  «  =  l4° ,  deflecting  the 

trailing  edge  flaps  15°  has  little  effect  on  ACp  ,  as  shown  in  Figures  19 
and  2b. 

Some  indication  of  the  effect  of  deflecting  the  leading  edge  flaps  8 
degrees  is  obtained  by  comparing  Figures  22  and  25.  The  greatest  effect  is 
to  raise  ACprms  at  the  aft  part  of  the  78  percent  semispan  by  a  factor  of 
2  to  3.  - 

4,2.3  Bending  Moment  Measurements  -  Wing  root  bending  moment  data  were 
analysed  to  help  assess  the  magnitude  of  the  buffet  dynamic  loads  problem. 

The  total  wing  root  bending  moment  signal  was  resolved  into  a  steady  compon¬ 
ent  and  an  unsteady  component.  The  rms  of  the  latter  signal  was  corrected 
by  subtracting  (on  a  mean  square  basis)  tunnel  background  noise  contributions 
yielding  the  rms  banding  moment  due  to  buffet.  These  rms  bending  moments 
%,  expressed  as  a. fraction  of  the  corresponding  steady  bending  moments, 

M,  are  presented  In  Figures  2o  through  20  for  Configuration  1  -  clean  wing. 
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(Examination  of  uncorrected  data  for  other  Mach  numbers  and  other  configura¬ 
tions  suggested  that  the  data  in  these  figures  encompasses  the  most  adverse 
a,,/M  ratios.) . 

The  maximum  rms  root  bending  moments  are  of  the  order  of  5-7%  of  the 
steady  moment,  being  maximum  at  about  MachO. 8.  Preliminary  unpublished  infor¬ 
mation  from  other  current  dynamic  buffet  investigations  yields  pretty  much 
the  same  story  -  maximum  ratio  of  cr^/M  of  about  the  same  order,  also  being 
most  critical  around  Mach 0.8.  The  orders  of  magnitude  coming  out  of  this 
and  other  recent  programs  suggest  something  about  the  magnitude  of  the 
buffet  dynamic  loads  problems .  Translated  into  operating  envelopes  of  fighter 
aircraft  (such  as  the  F-4),  cr^  translates '  roughly  as  5-7%  of  the  static  design 
load  at  sea  level;  but  under  typical  transonic  combat  maneuvering  at  30,000 
feet  for  existing  aircraft,  is  about  2-3%  of  the  design  static  load  at  the 
same  time  the  actual  static  load  is  roughly  only  1/3  to  1/2  of  design  static 
limit  load.  The  3  dynamic  dispersion  applied  to  the  static  load  is  then 
6-9%  of  the  static  load.  It  begins  to  appear  that  the  buffet  dynamic  loads 
problem  is  not  of  major  impact  on  static  design  criteria. 

4.3  Analysis  of  Spectral  Data 

Spectral  data  were  obtained  by  the  procedures  explained  in  Section  3 
and  data  are  presented  as  plots  of  power  spectral  density  and  normalized  cross 
power  spectral  density  versus  frequency.  The  buffet  test  conditions  chosen 
for  data  analysis  and  the  reasons  for  choosing  them  are  as  stated  at  the 
beginning  of  Section  4.2,2, 

4.3.1  Noise  Power  Spectral  Densities  -  The  noise  data  measured  at  CL  =  0 
and  discussed  in  Paragraph  4.2  were  also  analyzed  for  spectral  content.  Mach 
0.7  noise  data  at  four  points  along  the  70  percent  chord  line  are  shown  in 

Figures  29  through  32.  These  and  other  noise  data  indicate  a  strong  peak  in 

the  neighborhood  of  400  Hz  on  all  parts  of  the  wing  at  M  =  0.7.  In  most  of 
the  figures,  the  value  of  this  peak  is  between  2  x  10  '  and  4  x  10  (psi)  /Hz. 
Peaks  at  60,  120,  240  Hz,  etc.,  are  assumed  to  be  due  to  electrical  noise. 

Pressure  spectra  for  four  points  along  the  70  percent  chord  line  at 
Cji  -  0  and  M  =  0.9  are  shown  in  Figures  33  through  36.  These  and  other 

spectra  exhibit  fairly  broad  peaks  between  500  and  700  Hz  at  levels  generally 

between  2  x  10  ^  and  1  x  10  ^  (psi J4" /Hz. 


Two  spectral  plots  for  =  0  and'M  =  1.0-  are  shown  in  Figures  37  and 

38.  These  and  some  other  spectral  data  for  this  test  condition  exhibit  small  , 

broad  peaks  between  500  and  1000  Hz.  .The  levels  vary,  but  generally  they  - 
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are  below  10  (psi)  /Hz. 

*+.3.2  Buffet  Power  Spectral  Densities  -  Figures  39  through  6l  contain 
the  power  spectral  density  functions  for  the  fluctuating  pressures  at  Mach 
0.7  and  =  12°.  They  show  the  spatial  variation  of> the  fluctuating  pres¬ 
sure  spectra  over  the  wing  planform.  There  are  four  significant  character¬ 
istics  evident  to  varying  degree  in  these  plots:  a  broad  peak  centered  about 
650  hz ,  a  broad  peak  centered  about  1500  Hz ,  an  extremely  broad  peak  with  a 
maximum  at  650  Hz  with  delayed  roll-off  beyond  the  peak,  and  a  negative  slope 
from  about  10  to  200  Hz.  Transducers  whose  PSD  plots  exhibit  any  of  the  first 
three  characteristics  are  indicated  on  the  wing  planform  in  Figure  62.  Peaks 
at  650  Hz  are  evident  for  the  transducers  in  a  cone-shaped  region  emanating 
from  the  snag.  Based  on  flow  visualization  (oil  and  tuft)  studies  and 
analysis  of  the  static  pressures  over  the  wing,  this  general  region  is 
thought  to  be  affected  by  a  vortex  caused  by  the  snag.  Thus  the  650  Hz  peaks 
are  also  thought  to  be  caused  by  this  vortex..  The  transducers  along  the  78 
percent  span  line  have  particularly  broad  peaks  at  650  Hz  with  delayed  roll¬ 
off  beyond  the  peak.  The  spectrums  for  transducers  12  and  19  along  the  68 
percent  semispan  line  have  no  evidence  of  a  peak  at  650  Hz  although  the 
spectra  to  either  side  of  them  do  exhibit  the  peak. 

Transducers  whose  PSD's  have  1500  Hz  peaks  are  in  two  regions  as  shown 
in  Figure  62.  For  one  group,  it  is  thought  that  the  peak  is  due  to  the  lead¬ 
ing  edge  vortex.  For  the  other  group  the  peak  is  thought  to  be  due  to  the 
tip  vortex.  The  existence  of  these  vortices  is  discussed  in  Appendix  A. 
Estimated  boundaries  for  the  regions  affected  by  these  vortices  are  shown 
in  Figure  62.  The  estimates  are  based  on  the  positions  of  the  transducers 
whose  PBD's  exhibit  the  1500  Hz  peaks . 

The  frequencies  associated  with  the  different  voi’tices  are  thought  to 
be  related  to  some  kind  of  circulation  or  rotation  frequency  of  the  vortices. 

Negative  slopes  from  about  10  to  200  Hz  are  evident  in  the  PSD's  for  the 
transducers  along  the  tip  and  leading  edge.  Because  of  the  locations  of  these 
transducers,  it  is  likely  that  the  negative  slopes  are  also  related  to  the 
leading  edge  and  tip  vortices  although  in  a  different  manner  from  the  trans¬ 
ducers  with  i  500  Hz  peaks .  The  two  phenomena  do  not  always  occur  together. 


The  power  spectral  density  'curves 'in  Figures  39  through  6l  show  no 
indication- of  the  model  wing  -resonances  listed  in  Table  III..  The  only  model 
resonances  which  are  evident  in  these  curves  are  pitch  and  roll  of  the  model 
vibrating  on  the  balance.  Many  of  the  power  spectral  density  curves  contain 
small  peaks  at  the  frequencies  (12  and  18  Hz)  of  these  resonances  (Table  IV).- 

Power  spectral  density  plots  for  two  fluctuating  pressures  along  the  JO 
percent  chord  line  at  Mach  0.9  and  a =  8°  are  shown  in  Figures  63  and  64. 

These  plots  are  representative  ones  for  this  test  condition.  Most  of  the 
curves  have  negative  slopes  from  about  10  Hz  out  to  around  200  Hz.  The 
spectral  density  for  transducer  10  (Figure  63)  contains  a  low  peak  at  1500  Hz 
and  for  transducer  11  (Figure  64)  there  is  a  medium  peak  around  1000  Hz. 

These  frequencies  and  those  exhibited  by  the  power  spectral  density  plots  of 
the  other  pressure  transducers  are  indicated  on  the  wing  planform  in  Figure  65 
While  the  situation  here  is  not  as  clear  as  it  was  at  Mach  0.7  and  aw  =  12°, 
regions  affected  by  three  vortices  are  evident  and  are  indicated  in  Figure  65. 
An  inner  wing  leading  wing  vortex  is  associated  with  peaks  in  the  neighborhood 
of  1500  Hz  at  transducers  8,  12,  19,  and  20.  The  snag  vortex  is  associated 
with  1000  to  1200  Hz  peaks  measured  by  chordwise  transducers  along  the  78  per¬ 
cent  semispan  line.  Chordwise  transducers  along  the  93  percent  semispan  line 
have  peaks  from  1200  to  1500  Hz.  There  is  some  question  as  to  what  vortex 
causes  them.  If  the  peaks  were  not  evident  for  transducer  1,  it  would  be 
assumed  they  are  associated  with  the  tip  vortex.  Since  a  peak  is  evident  for 
transducer  1,  it  is  indicated  in  Figure  65  that  the  peaks  are  due  to  an  outer 
wing  leading  edge  vortex.  Since  this  is  a  unique  case  among  the  test  condi¬ 
tions  analyzed,  this  vortex  is  included  with  the  tip  vortex  in  the  prediction 
method . 

Power  spectral  density  curves  for  four  fluctuating  pressures  measured 
with  transducers  located  along  the  JO  percent  chord  at  Mach  0.9  and  av  =  12° 
are  shown  in  Figures  66  through  69.  The  curves  for  transducers  10  and  11  are 
shown  in  Figures  66  and  67 »  respectively »  and  exhibit  strong  peaks  at  425  Hz. 
The  curves  for  transducers  13  and  l4  are  shown  in  Figures  68  and  69.  These 
curves  have  peaks  at  750  Hz.  A  summary  of  the  peaks  in  the  pressure  power 
spectral  densities  at  this  test  condition  is  presented  in  Figure  70,  and  the 
areas  of  the  wing  over  which  the  vortex  frequencies  are  evident  in  the 
fluctuating  pressures  are  indicated, 
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y  Figures  71  through  73  contain  power  spectral  density  curves  for  fluctu-. 
ating  pressures  measured  with  transducers  10,  11,; and  13  at  Mach  0.9' and 
The  curve  for  transducer  10  does  not- exhibit  a  significant  high 
frequency,  peak.  The; curve  for  transducer'll  'has  peaks  at  375  and  700  Hz. 

For  transducer;. 13,  there  is  a  peak  only  at  J00  Hz.  As  shown  in  Figure  74,  ’ 
.the  375  Hz  frequency  is-  associated  with  the  snag  vortex,  and  the  J00  Hz 
frequency  is  associated  with  . the  leading  edge  and  tip- vortices . 

y.  Power, spectral .density  plots  for  four  fluctuating  pressure  transducers  - 
along  the  70  percent  chord  at  Mach  l'and  aw  =  12°  are  shown  in  Figures  75 
through  78." -  .No  trend  is  apparent  in  these  and  the  other  plots  examined  at 
-.this  test  condition.  A  possible  explanation,  for  the  lack  of.  a  pattern  is 
that-  -the  spectrums  may  be  affected  by  shock- waves  on  the  wing. 

■"The  general- trends  discussed  above  for  configuration  1  in  different  flow 
conditions  are  also  apparent  for  configurations  5  and  6.  Some  PSD  plots  for 
-  fluctuating  pressures  measured  along  the  70  percent  chord  line  and  drawings 
showing  how  peaks  in  these  and  other  plots  are  located  on  the  wing  are  shown 
in  Figures  79  through  96 . 

it,  3. 3  Hatio  of  Peak  to  Low  Frequency  Spectral  Levels ,  -  The  model  for 
the  spectral  shape  of  the  fluctuating  pressures  developed  in  Section  5  has  as 
one  of  its  parameters  a  quantity  r\ ,  which  is  the  ratio  of  the  magnitude  of 
the  peak  level  in  the  spectrum  to  the  low  frequency  level.  The  low  frequency 
level  is  assumed  constant  from  a  very  low  frequency  out  to  the  frequency  at 
which  the  spectrum  starts  to  increase  towards  the  peak.  Values  of  r>  have 
been  computed  for  the  buffeting  pressure  power  spectral  density  plots  which 
contain  significant  peaks  associated  with  the  vortices  discussed  in  Section 
4,3,2,  These  values  are  plotted  as  functions  of  percent  chord  in  Figures  97 
through  104.  For  plots  containing  more  than  one  such  peak,  n  was  computed 
for  the  higher  one.  For  plots  which  are  not  constant  at  the  low  frequencies, 
an  average  value  of  the  power  spectral  density  at  the  low  frequencies  was 
used  in  the  denominator  of  n ,  When  this  average  value  is  greater  than  the 
value  of  the  peak,  then  U  <1. 

There  are  two  patterns  evident  in  the  plots  of  n  versus  percent  chord  in 
Figures  97  through  104.  The  first  is  that  n  is  generally  larger  along  the 
78  percent  seraispan  line  than  along  the  other  constant  span  lines  on  which 
measurements  were  token.  This  pattern  was  also  evident  in  the  plots  of 
overall  rms  pressures  in  Figures  l4  through  25*  This  semispan  line  is  close 


.to  the  center  of  the  region  affected  by  the  snag  vortex.  The  second  pattern, 
is  that  n  tends  to  increase  from,  the  leading  edge  to  the  trailing  edge,  of 
the  wing. 

Buffet  Cross  Power  Spectral  Densities  -  Normalized  cross  power 
spectral  densities  were  computed  to  determine  (a)  whether  a  convection  veloc¬ 
ity  exists  in  the  separated  region  (and  if  so  its  value),  and  (b)  the  value 
of  the  spatial  decay.  The  functions  were  plotted  as  magnitude,  y,  and,  phase, 
0,  versus  frequency.  The  spatial  decay  parameter  is  determined  from  the 
value  of  the  magnitude  as  a  function  of  the  distance  between  the  transducers 
in  the  pair.  If  a  condition  of  pure  convection  independent  of  frequency 
exists,  then  the  phase  is  proportional  to  frequency.  The  normalized  cross 
power  spectral  density  plots  in  Figures  105  through  12lt  are  samples  of  many 
which  show  that  this  is  not  the  case  for  the  F-h  wing  as  a  whole  in  separated 
flow.  These  plots  are  for  pressures  measured  with  transducer  pairs  (6,  10) 
and  (12,  18)  and  the  test  conditions  for  configuration  1  covered  in 
Section  h.3.2  in  the  discussion  of  power  spectral  densities. 

Figures  105  and  106  show  the  magnitude  and  phase,  respectively,  of  the 
normalized  cross  power  spectral  density  function  for  transducer  pair  (6,  10) 
at  Mach  0.7  and  =  12°.  There  are  large  changes  in  y  between  700  and  1000 
Hz  and  between  1300  and  1700  Hz,  An  average  value  for  y,  taking  into 
account  the  logarithmic  frequency  scale,  is  about  0.35.  This  and  similar 
values  are  used  to  determine  the  spatial  decay.  The  phase  has  varying  slopes, 
but  the  trend  is  positive  from  0  to  around  900  Hz.  From  900  Hz  to  2000  Hz, 
the  trend  in  slope  is  negative.  Thus,  the  phase  is  definitely  not  proportiona 
to  frequency  in  this  case.  Figures  107  and  108  present  y  and  0  for  pair 
(12,  18).  Again,  y  is  not  constant;  an  approximate  average  value  is  0.35. 

The  phase  again  is  not  proportional  to  frequency,  although  it  is  more  nearly 
so  than  that  for  pair  (6,  10). 

Figures  109  and  110  contain  the  y  and  0  functions  for  pair  (6,  10)  at 
Mach  0.9  and  a  »  8°.  The  relatively  small  variation  in  y  over  most  of  the 
frequency  band  and  the  good  approximation  of  0  to  a  linear  function. of 
frequency  indicate  the  existence  of  a  nearly  pure  (i.e.,  independent  of 
frequency)  time  delay  from  transducer  10  to  transducer  6.  One  explanation 
for  the  time  delay  is  the  existence  of  a  local  convection  velocity  having 
its  fore-aft  component  directed  forward.  The  time  delay  might  also  bo 
partially  caused  by  a  difference  in  propagation  distances  from  a  fluctuating 


pressure  source  above  the  wing.  Figures  ill  and  112  show  y  and  6  for  pair 

(12,  18) .  There  is  a  large  variation  in  y.  The  graph  of  phase  versus 

frequency  has  three  distinct  average  slopes.  Therefore,  there  is  no  pure 

time  delay  between  transducers  11  and  18  and  no  pure  convection. 

The  plots  of  magnitude  and  phase  of  the  normalized  cross  power  spectral 

density  for  pair  (6,  10)  at  Mach  0.9  and  a  =  12°  are  shown  in  Figures  113 

w 

and  lilt  and  indicate  the  absence  of  a  pure  time  delay.  The  plots  for  pair 
(12,  18)  are  shown  in  Figures  115  and  ll6.  They  also  indicate  the  absence 
of  a  pure  time  delay.  The  snag  vortex  frequency  of  425  Hz  is  evident  in 
Figure  115. 

Plots  of  y  and  8  for  pairs  (6,  10)  and  (12,  18)  at  Mach  0.9  and  =  l4° 

are  shown  in  Figures  117  through  120.  Figure  118  again  shows  that  there  is 

no  pure  time  delay  between  transducers  6  and  10.  The  phase  plot  for  pair 

(12,  18)  has  several  regions  of  different  average  slopes. 

The  normalized  cross  power  spectral  densities  for  pairs  (6,  .10)  and 

(12,  18)  at  Mach  1  and  a  =  12°  are  shown  in  Figures  121  through  124.  The 

w 

plots  of  y  show  that  the  signals  in  the  first  pair  are  more  strongly  corre¬ 
lated  than  are  those  in  the  second  pair.  The  plot  of  0  for  pair  (6,  10)  has 
a  change  in  average  slope  at  about  1000  Hzj  this  again  shows  the  absence  of 
a  pure  time  delay  between  transducers  6  and  10,  The  plot  of  0  for  pair 
(12,  l8)  also  shows  the  absence  of  a  pure  time  delay. 

Examination  of  many  cross  power  spectral  density  functions ,  of  which 
the  ones  discussed  above  are  samples,  has  led  to  the  conclusion  that  con¬ 
vection  is  not  a  highly  significant  parameter  in  buffet  of  highly  swept 
fighter  wings  and  that  convection  effects  can  be  neglected  in  the  development 
of  a  first  generation  mathematical  model  for  the  prediction  of  buffet  loads 
The  observations  which  led  to  this  conclusion  were  the  following: 

o  The  correlation  between  the  fluctuating  pressures  at  two  points 
decays  fairly  rapidly  with  distance  between  the  points, 
o  The  existence  of  convection  is  indicated  when  the  phase  el*  the  cross 
power  spectral  density  function  is  linear  with  frequency  ?uid  passes 
through  the  origin.  In  general ,  the  data  did  not  follow  this  trend. 
o  The  interaction  of  vortices  originating  from  the  leading  edge,  snag, 
ami  tip  is  80  complex  that  convection  in  a  broad  sense  seems  to  be 
fairly  random  in  nature.  Convection  appears  to  be  highly  localised. 


Additional  justification  was  provided  by  cursory  a  study  made  to  establish 

the  general  effects  of  n  localized  but  independent  convection  zones.  The 

results  indicated  that  the  structural  response  of  a  simple  torsion  model 

would  be  reduced  by  a  factor  ^  compared  to  the  response  for  a  pure  convection 

situation.  The  flow  visualization  oil  and  tuft  studies  indicate  that  n  is  6 

or  greater;  hence,  even  if  the  localized  convection  effects  were  orderly, 

little  convection  influence  on  response  would  be  expected. 

k,3.5  Tail  Buffet  -  For  an  examination  of  tail  buffet,  power  spectral 

densities  were  calculated  for  transducer  27,  located  on  the  left  horizontal 

tail  (stabilator)  at  Mach  0.9  and  aw  =  8°,  12°,  and  l4°.  These  curves  are 

shown  in  Figures  125  through  127.  The  corresponding  curves  for  transducer 

21,  located  on  the  left  trailing  edge  flap,  are  shown  in  Figures  128  through 

130  for  comparison.  At  =  8°,  the  fluctuating  pressures  on  the  stabilator  . 

were  due  to  tunnel  noise.  At  &w  =  12°,  the  power  spectral  density  of  the  , 

fluctuating  pressures  there  has  a  peak  at  the  snag  vortex  frequency,  but  the 

leading  edge  vortex  frequency  is  not  evident.  The  overall  rms  level  on  the 

stabilator  is  lower  than  that  at  most  points  on  the  wing.  Observations  for 

=  12°  apply  equally  well  for  ay  =  lh°. 

The  normalized  eross  power  spectral  density  between  transducer  21  on  the 

trailing  edge  flap  and  transducer  27  on  the  stabilator  at  Mach  0.9  and  etw  »  12° 

is  shown  in  Figures  131  and  132.  This  and  similar  functions  at  Mach  0.9  and 

u  »  8°  and  ik  show  that  there  is  little  correlation  between  the  fluctuating 
w 

pressures  on  the  wing  and  those  on  the  tail.  There  is  some  coherence  at  the 

snug  vortex  frequency  at  u  =  12°  said  it0 .  This  indicates  the  existence  of 

some  kind  of  convection  from  the  wing  to  the  tail,  but  not  one  with  which 

there  is  associated  a  clear  convection  velocity. 

The  erois  power  spectral  density  between  the  pressures  on  either  aide 

of  the  stabilator  at  Mach  0.9  and  a,,  =  12°  was  computed.  As  was  expected, 

w 

there  is  essentially  no  correlation  between  these  pressures. 
it.lt  Comparison  of  Wind  Tunnel  and  Flight  Data 

The  comparison  of  wind  tunnel  (model)  and  flight  (full  scale)  fluctuating 
pressures  in  studied  in  two  parts;  (l)  the  overall  rms  level,  and  (2)  the 
•spectral-  shapes.  These  parts  are  considered  separately  using  commonly 
accepted  relationships  for  fluctuating  pressures  (tic*?reneea  18,  19,  20). 

Flight  data  are  from  a  VF-4E  airplane  flight  teat -program  which  was  reported 
If.  He fo reuse  2. 


b.b.l  Overall  rms  -  The  overall  rms  level,  a,  of  fluctuating  pressures 
is  generally  assumed  to  be  proportional  to  the  freestream  dynamic  pressure,  q. 
Hence,  the  rms  pressure  coefficient,  A  On  ~  c/q,  is  considered  to  be  inde- 
pendent  of  scale.  Comparisons  of  model  and  full  scale  rms  pressure  coefficient 
variations  with  angle  of  attack  at  Mach  0,8  and  0.9  are  presented  in  Figures 
13  3  through  lb 0.  Full  scale  data  are  show.:  for  flights  at  25,000  feet  and 
35,OOG  feet.  The  comparisons  are  for  transducers  17  (—  =  0.9,  =  0.86), 

18  (f  *  0.9,  «  0.78) ,  21  (-  =  0.9,  =  0.33) ,  and  2b  (|  =  0.95, 

=  0.33).  At  Mach  0.8,  the  model  and  full  scale  curves  for  transducer  17 
(Figure  13  3)  have  different  curvatures  and  compare  unfavorably.  For  trans¬ 
ducer  18  (Figure  13  b),  the  curves  are  very  similar  in  shape  and  the  curve  for 
flight  data  is  displaced  2°  to  3°  to  the  right  of  that  for  the  tunnel  data. 

A  given  value  of  ACp^_,  would  therefore  occur  at  an  angle  of  attack  2°  to  3° 
lower  in  the  tunnel  than  in  flight.  For  transducer  21  (Figure  135)  the 
curves  for  model  and  full  seale  are  similar.  To  compare  values  of  A^’prrn3  at 
a  given  angle  of  attack,  the  corresponding  noise  values  (those  fit  low  angle 
of  attack)  should  be  subtracted  on  a  mean  square  basis.  When  Acprms  is  large 
compared  to  its  corresponding  noise  level,  however,  the  effect  of  the  noise 
is  negligible.  Thus,  at  the  higher  angles  of  attack  shown,  the  model  data  *s 
not  only  higher  than  the  full  scale  data  but  it  is  also  increasing  faster. 
Transducer  2b  is  near  transducer  21  and  the  data  (although  limited  for  flight) 
for  transducer  2b  (Figure  136)  appear  essentially  the  same  as  those  for  transducer  21. 

At  Mach  0.9,  the  model  and  full  scale  curves  for  transducer  17  (Figure  137) 

are  similar;  however,  above  5®  the  model  data  are  higher  and  above  12°  they 

are  increasing  faster  than  are  the  full  scale  data.  For  transducer  18,  the 

curves  for  tunnel  data  are  much  higher,  particularly  from  uw  =  8°  through  lb°. 

Figure  139 again  shows  that  the  flight  ana  tunnel  curves  ere  similar.  The 

differences  are  that  the  curve  for  the  tunnel  in  higher  and  that  from  ll0  to 

the  curvatures  of  the  two  curves  are  different.  The  limited  amount  of 

flight  data  for  transducer  2b  at  Mach  0,9  agrees  fairly  well  with  the  tunnel 

data  for  that  transducer  except  that  the  tunnel  data  are  lower.  This  is  shown 

4  **  _  ■ 

■  in  FigurelbO,  '  ' 

in  summary,  the.ACp^  versus  <*w  curves  for  model  and  full  scale,  which 
by  eecseon  assumption  should  coincide,  ai*e  similar  but  are  not  identical. 

%r  transducers  21  and  2b,  located  on  the  trailing  edge  flu-  ,  buffet  AC™ 
is  higher  in  the  tunnel  than  in  flight.  It  is  thought  that  some  of  the 
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differences  between  the  tunnel  and  flight  curves  for  transducers  IT  and  18 
can  be  explained  by  the  strong  pressure  field  generated  by  the  snag  vortex. 

If  the  axis  of  this  vortex  (and,  hence,  other  parts  of  it)  in  in  a  slightly 
different  position  relative  to  these  transducers  in  the  wind  tunnel  than  it 
is  in  flight,  then  ACp_^s  could  vary  considerably  between  tunnel  and  flight. 

The  strongest  evidence  for  such  a  possibility  is  shown  in  Figure  134.  A 
given  value  of  is  reached  in  the  tunnel  at  a  2°  to  3°  lower  angle  of 

attack  than  in  flight.  If  the  axis  of  the  snag  vortex  reaches  a  given 
position  relative  to  transducer  18  at  a  2°  to  3°  lower  angle  of  attack  in  the 
tunnel  than  it  does  in  flight ,  then  Figure  134  would  be  explained. 

Another  possible  explanation  for  some  of  the  differences  between  the 

ACL  versus  a  curves  for  wind  tunnel  and  flight  is  an  error  in  the  flight 

measurements  of  due  to  nose  boom  upwash.  If  such  an  error  does  exist,  the 

values  of  a  (flight)  shown  are  thought  to  be  high  by  about  1  degree, 
w 

4.4.2  Frequency  Distribution  -  The  key  element  in  comparing  the 

spectral  shapes  of  the  fluctuating  pressures  in  the  tunnel  and  in  flight  is 

the  commonly  accepted  assumption  that  the  frequency,  f,  associated  with  a 

pressure  phenomenon  is  proportional  to  the  freestream  velocity,  V,  and  inversely 

proportional  to  a  reference  length,  b.  In  other  words,  the  non-dimensional 
flo 

frequency,  k  =  —  ,  is  assumed  independent  of  scale.  One  problem  is  choosing 
a  suitable  reference  length.  In  boundary  layer  studies,  the  reference  length 
is  often  chosen  as  the  displacement  thickness  or  the  momentum  thickness,  In 
the  separated  flow  associated  with  buffet,  there  is  no  classical  boundary 
layer  as  such.  It  will  be  shown  that  a  reference  length  associated  with  the 
structure  (wing,  for  example)  gives  good  scaling  results. 

For  the  comparison  of  frequency  distributions  -  i.e.,  of  nondimen3ional 
spectral  shapes  rather  than  of  actual  PUD  levels  -  flight  data  were  scaled 
to  model  data  using  the  following  equations: 


^F8 

f1  a  ■■■■  -  ».m  — 4* 

K  VM  M 

n  9  b  V 

«F3  FS  M 


(4-1) 


(4-2) 


whore 


0  is  power  spectral  density  of  fluctuating  pressure,  (psir/Ilz 
f  is  frequency,  Ha 
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o'  is  overall  rms  fluctuating  pressure 
b  is  structural  reference  length 
V  is  freestream  velocity 
M  (subscript)  denotes  model 
FS  (subscript)  denotes  full  scale 

Equation  (4-1)  assumes  that  ~  is  independent  of  scale.  Equation  (4-2)  follows 
from  equation  ( 4—1 )  and  the  fact  that  a  =  /”  G(f)  df. 

Wind  tunnel  and  flight  pressure  spectral  shapes  are  compared  in  Figures  lit 
through  146.  The  flight  spectral  shape  was  obtained  by  applying  equations 
(4-1)  and  (4-2)  to  power  spectral  density  data  which  were  plotted  on  a  linear- 
linear  graph.  A  scries  of  points  near  zero  on  the  linear  graph  appear  on  the 
log-log  graph  as  a  straight  line.  The  power  spectral  density  functions  for 
the  flight  data  were  computed  digitally  at  2  Hz  intervals.  The  equivalent 
bandwidths  and  averaging  times  were  4  Hz  and  2.5  seconds,  respectively.  Some 
points  from  the  linear-linear  plots  were  not  plotted  in  Figures  l4l  through 
146  because  they  were  not  needed  to  display  the  principal,  characteristics 
of  the  spectral  shapes . 

The  six  comparisons  of  model  and  full  scale  spectral  shapes  shown  in 

Figures  ]4l  through  146  are  for  transducers  17  (~  ~  0.90,  =  0.36)  and 

l8  (—  =  0.90,  ~  0.78)  and  three  test  conditions.  The  test  conditions  are 

c  b 

shown  in  Table  IX.  All  sir  comparisons  show  good  agreement,  particularly 


Table  IX 

Test  Conditions  for  Comparison  of  Model  and 
Full  Seale  Spectral  Shapes 


Configuration 

Much  No. 

Wing  Anglo  of  Attack,  Dog 

1 

0,8 

10 

1 

0.80 

10 

0 

0.8 

0 

in  the  regions  of  the  peaks  mured  by  vortlct?*;,  Pom*  of*  the  flight  spectral 
shapes  contained  abroad  peak  around  90  (sealed  to  model  frequency)  that 
not  evident  in  the  tunnel  spectral  shape.  Thin  In  in  the  neighborhood  oi 
the  (scaled)  first  hu tiding  frequency  of  the  full-scale  wing  find,  coinciden- 
tnliy,  of  the  first  bending  frequency  of  the  model  wing  alee.  The  fact  tha* 
the  broad  high  frequency  peaks  for  the  model  and  scaled  full  scale  spectral 
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SECTION  5  ■?: 

PREDICTION  METHOD  FOR  FLUCTUATING  PRESSURES 

i 

The  buffet  prediction., method  described  herein  requires  that  spatial, 
magnitude  and  spectral  characteristics  of  random,  fluctuating  pressures  be  • 
described  on  the  upper  wing  surface,  where  significant  flow  separation 
occurs  at  the  higher  angles  of  attack  in  the  buffet  regime.  This  section 
presents  a  general  discussion  of  the  nature  of  the  flow  pattern  over  a  highly 
swept  wing  in  buffet,  the  nature  of  the  fluctuating  pressures  and  their 
distribution  on  the  wing,  and  a  feasible  method  for  prediction.  A  mathemati¬ 
cal  model  is  developed  and  a  step-by-step  method  is  set  forth  for  using  the 
pressure  prediction  method.  ■ 

5.1  General 

It  was  shown  in  Sect, ion  h  and  in  Appendix  A  that  the  aerodynamic  flow 
over  the  F-h  wing  is  extremely  complex,  and  that  the  associated  fluctuating 
pressures  are  similarly  complex.  'Two  factors  have  dominated  the  study  of 
the  nature  of  fluctuating  pressures  to  be  applied  to  an  upper  wing  surface: 

(1)  Some  rms  pressure  levels  and  spectral  peaks  have  been  related  to 
vortices  (Section  and  Appendix  A)  but  a  comprehensive  pattern 
describing  the  flow  on  the  F-4  wing  has  not  emerged, 

(2)  In  the  preliminary  design  stage  of  a  new  airplane,  knowledge  of 

wing  flow  characteristics  in  buffet  is  extremely  limited  because 
the  available  tools  (lift  and  drag,  data,  flow  visualization  data) 
have  not  been  sufficiently  developed  to  provide  reliable  correla¬ 
tions.  •  :  . 

Therefore,  a  significant  portion  of  the  study  effort  was  devoted  to  examining 
the  data  obtained  with  F~Jt  wing  studies  for  the  purpose  of  establishing 
generalisations  that  are  believed  to  be  common  to  all.  wings. 

It  was  originally  enviaiov-d  that  the  mathematical  modeling  of  the 
input  buffeting  pressures  eou.u  be  based  on •  .spectral  functions  formulated  in 
terms  of  ms  value,  spectral  shape,  length  scale  value  L,  and  convection 
velocity  Vf>.  The  experimental  data1  do  not  indicate  an  overall  convection 
pattern  m  the  F-li  wings ;  if  "convective!*  phenomena  are  present,  they 
probably  occur  only  in  highly  localized  regions.  Attempts  to  extract  V 
and  L  have  little  meaning,  therefore,  and  characterization  in  terms;  of  more 
general  spectral  and  cross-smctrul  functions  was  adopted.  It  is  significant 


to  understand  that  the  direct  use  of  the  experimentally  indicated' cross- 
spectral  functions  does  not  exclude  convection  effects  or  imply  that  some  • 
important  factor  is  being  overlooked.  Convective  effects ,  if  -they  are 
present,  are  implicit  in  these/  functions.  The-  practical  use. of  spectral'--:': 
functions  always  involves  the  ratio  L/V:^;  thus  in  general,  it  is  .-only  necessary 
to  find  the  value  of  this  ratio.  The  separation  into  values  L  and  V  •  is."’7  v 
meaningful,  and  serves  as  a  convenient  way  to  describe  the  flow  in  a-. more 
detailed  fashion,  only  when  convection  is  clearly  .present, 

5.2  Nature  of  Pressure  Function 

The  characteristics  of  the  buffeting  fluctuating  pressures  that  are  of 
importance  in  computing  structural  response  are  the  -power  spectra  of  the,, 
pressures  at  various  locations  on  the  wing  surface,  and  the  cross,  spectra 
between  separated  points.  It- -is  convenient  to  split  the  representation  of 
each  of  these  functions  into  two  parts,  the  rms  values  of  the  pressures, 
and  the  spectral  function  shape.  The  variation  of  these  quantities  with 
angle  of  attack  a  and  Mach  number  M  forms  the  essence  of  the  method. 

5.2.1  rms  Levels  -  Generally,  the  descriptions  of  the  rms  pressure 
fluctuating  levels  in  shear  flow,  turbulent  boundary  layers,  or  buffeting 
is  commonly  and  conveniently  given  in  terms  of  a  reference  q.  This  conven- 
tion  has  also  been  used  here  and  the  equation  for  overall  rms  pressure 
level  that  has  been  adopted  is 


where 


q  is  freestream  dynamic  pressure 


(5-1) 


AG  is  the  local  fluctuating  pressure  coefficient. 

^rms 

Application  of  this  equation  requires  that  the  coefficient  ACPmg  be  known 
at  various  locations  on  the  wing  and  for  various  values  of  and  M. 

Simplified  models  representing  the  variation  of  ACp^  over  the  wing 
surface  were  derived  by  considering  the  spatial  distribution  of  the  data 
contained  in  Figures  Ik  through  25.  Aside  from  this  data  itself,  several 
other  factors  were  considered  in  arriving  at  these  representations,  such  os 
those  listed  below: 


'The  applicability  of  the  prediction  method  would  be  limited  to  the 

:  ,7^;  4etermi nation  -of- ^he,  primary  loads  in  the  structure  in  contrast  to 

-  localized  loads  such  as  on  individual  wing  skin  panels.  Also  the 

corresponding  primary  structure  model  deformations  influencing  these 
-loads  do  not  vary  greatly  over  localized  areas.  Hence,  detailed 
definition  of  localized  fluctuating  pressure  variations  is  an 
unnecessary  luxury. 

(b)  Magnitudes  of  fluctuating  pressures  and  structural  loads  were 

a  fraction  of  the  corresponding  static  values.  Hence,  any  rea¬ 
sonable  rough  approximation  of  the  fluctuating  pressures  is 
consistent  with  its  overall  structural  impact. 

(c)  The  intent  of  the  research  program  was  to  develop  a  first  genera¬ 
tion  method  for  buffet  load  prediction  during  preliminary  design. 
Excessive  detail  and  complexity  would  be  antithesis  to  this  intent 
and  could,  in  fact:-,  discourage  its  use  for  that  purpose. 

(d)  In  spite  of  the  overall  complexity  of  the  flow  details,  a  gross 
behavior  pattern  seemed  to  emerge.  Excluding  the  singular  snag 
effect,  flow  separation  and  fluctuating  pressure  intensities 
tended  to  increase  outboard  along  the  span  and  rearward  along  the 
chord, 

(e)  The  present,  sbate-of-the-art  in  predicting  buffeting  pressures  on 
lifting  surfaces  is  yet  in  its  infancy.  The  limited  information 
available  is  inconclusive  beyond  first  order  effects. 

In  light  of  the  above  considerations,  it  was  decided  to  employ  a  linear 
variation  of  ^^pms  over  the  wing  such  as  those  shown  in  Figure  1U7.  An 
example  comparison  between  the  linear  model  of  Figure  ll*7(a)  and  measured 
data  points  from  Figures  15*  18,  and  19  (for  the  clean  wing)  is  shown  in 
Figure  lU£.  The  models  in  Figure  li*7  are  for  use  in  computational  response 
studies  and  constitute  the  ACpiina  models  to  be  used  for  wings  in  general. 

For  design  purposes  one  would  be  justified  in  designing  to  the  most  adverse 
environment.  Jn  linearized  representation.  Figure  lVf(a)  con  be  used  to 
this  f-t.  i  for  cl-  an  (ur;'"  upped)  configurations  while  Figure  ll»7(c)  can  be 
used  for  deflected  flap  configurations ,  at  least  for  the  present  state  of 
knowledge  in  this  area.  These  are  consistent  with  the  more  intense  environ¬ 
ments  encompassed  in  the  wind  tunnel  test  program  -  generally  the  high 
sue  sonic  Mach  numbers  and  angles  of  attack  from  about  10°  to  14°.  At  sonic 


and  supersonic  speeds  a  marked  reduction  in  fluctuating  pressures  was 
clearly  evident  (compare  Figure  llfT(b)  with  lV[(a)  for  example).  There 
was  some  evidence  of  a  leveling-off  trend  in  buffet  intensity  at  the  higher 
angles  of  attack  around  10°  -  lh° ,  but  such  data  were  too  sparse  to  be 
really  conclusive. 

5.2.2  Spectral  Shape  -  In  general  the  power  spectra  of  the  buffeting 
pressures  measured  on  the  F-4  wings  indicated  either  a  flat  or  white  noise 
spectral  shape  or  a  shape  with  a  pronounced  peak  at  some  discrete  frequency. 
These  shapes ,  corroborated  by  the  fact  that  the  autocorrelation  functions 
generally  were  of  the  decaying  oscillatory  wave  type*  suggest  the  adoption 
of  the  following  mathematical  description  of  the  autocorrelation  function : 


R(  t)  =  q~  e"8^  T I 


cos  (j)  i 
o 


(5-2) 


This  function  leads  to  a  wide  latitude  in  spectral  shapes.  For  =  0,  a 
spectrum  is  obtained  which  is  flat  in  the  low  frequency  range;  for  a  =  0,  it 
yields  a  spike  type  spectrum  at  the  frequency  uq.  Various  combinations  of 
a  and  ioq  yield  spectral  shapes  with  a  peak. 

A  non-dimensional  form  of  equation  (5-2)  is 


P  1  I 

R(s)  =  a  e  m  ^  cos  27rkos 


(5-3) 


where  m  =  a  constant  (empirical) 


f  c 

k  = 
o  V 


s  =  ~  =  non-dimensional  time 


f  =  frequency  at  which  peak  occurs 


c  =  reference  chord 


V  =  free  stream  velocity 

The  coefficient  "a"  in  equation  (5-2)  relates  to  equation  (5-3)  by 


r.tU; 


'"  v" >•  •*•'••.  :’  ^  c';yr^»^  '  ;v .  s.S..^  . :-'i  (.  -$b.  (?  ':‘r V. ••’>'":  1 


Since  m  is  a  constant,  it  is  noted  that  "a"  in  effect  represents  the 
ratio  Vc/L  previously  discussed;  V -takes,  the  place  of  V  ,  and  me  takes  the 
place  of  L. 

The  power  spectral  function  indicated  by  equation  (5-3)  is 


G(k)  _  m _ 1 _ +  _ 1 

a2  7r.  1  +  m2  (k  +  k)2  1  +  m2  (k  -  k)2 


b-k) 


which  contains  the  variable  parameters  m  and  kQ  (see  Reference  16,  pages 
86  and  87).  It  may  also  be  expressed  in  non-dimensional  form  as 


where 


°  2  it  0  .2 

«  !  +  x2  (1  ♦  ) 
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1  +  X2  (1  -  |~)‘ 
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(5-5) 


X  =  ink 


G(k)  =  r  G(f) 


i 
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and  the  ratio  k/kQ  is  a  plotting  variable. 

Various  spectra  indicated  by  equation  (5-5)  are  shown  in  Figure  1U9. 
5.2,3  Fitting  Data  to  the  Model  -  The  characteristics  of  the  mathemati¬ 
cal  model  that  can  be  directly  related  to  similar  characteristics  of  an 
experimental  data  plot  include 

0  the  location  of  the  peak  amplitude  (k  =  kQ) 

0  the  ratio  of  peak  amplitude  (at  k  a  kQ)  to  the  amplitude  at  zero 
frequency  (k  =  0). 

The  latter  ratio  is  called  n  and  can  be  calculated  in  terms  of  X  from 
equation  (5-5).  For  the  peal;  at  k  =  kQ, 

cog  -  A ( — L-$+1) 

^  V  1  +  k  a 

and  for  k  «  0 

■  irr  <777) 

0  1  +  X 
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and 
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"-GloT-  1+Ux2  = 


Solving  for  X  yields 


X  =  |  [4n  -  3  +  (l6n2 3  -  i6n  +  l)1/2]1/2,  n  >  0.933  (5-6) 


A  plot  of.  X  vs.  n  obtained  from  equation  (5-6)  is  shown  in  Figure  150. 

Values  of  n  <  1.0  have  no  physical  significance  in  the  mathematical  model 
and  are  not  included  in  the  figure, 

A  method  for  simplifying  and  systemizing  the  procedure  for  fitting 
equation  (5-5)  to  experimental  data  is  described  in  the  following  procedure: 


(l)  Determine  f  at  the  peak,  then 


k  =  f  c/V 
o  o 

(2)  Find  the  ratio  n  of  the  peak  (fQ)  value  to  the  f  =  0  value. 

(3)  Enter  Figure  150  with  this  n  and  read  out  X. 

An  example  of  how  the  spectrum  so  obtained  fits  a  measured  spectrum  is 
shown  in  Figure  151,  which  was  obtained  for  transducer  11  at  Mach  0.7 

and  a  =  12°' (see  Figure  62).  The  constant  values  used  in  the  mathematical 

w  _ 

model  were  fQ  =  700  Hz,  c  =  1.60  feet,  V  =  776  ft/sec,  o  =  1.1  psi, 
n=  12.5,  and  X  =  4,85. 

Based  on  the  results  of  Section  4,  the  existence  of  peaks  in  the 
power  spectral  density  of  the  fluctuating  pressures  is  dependent  on  the 
vortices  existing  over  the  wing.  Figures  62,  65,  70,  74,  82,  87,  92,  and 
96  indicate  on  the  basis  of  spectral  content  that  there  are  three  vortices 
of  interest.  These  are  associated  with  the  inner  wing  leading  edge,  the 
tip,  and  the  snog  of  the  wing.  The  frequencies  associated  with  these  vortices 
and  shown  in  the  figures  given  above  are  plotted  versus  wing  angle  of 
attack,  Mach  number,  and  configuration  in  Figure  152.  The  plots  versus 
angle  of  attack  are  for  configurations  1  and  5  at  Mach  0.9.  They  show 
chat  for  all  three  vortices  the  frequency  decreases  as  wing  angle  of 
attack  increases  over  the  range  from  8°  through  l4°  for  configuation  1. 

The  plots  versus  Mach  number  show  that  for  configuration  1  the  peak 
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frequencies  decrease  as  Mach  number  increases  over  the  range  0.7  through 
0.9.  For  configuration  5,  they  show  that  the  peak  frequencies  are  less 
dependent  upon  Mach  number  over  this  range.  As  mentioned  in  Section  4.3.2, 
no  pattern  emerged  in  the  spatial  variation  of  the  power  spectral  density 
plots  for  the  Mach  1  test  'condition.  -No  relationship  between  any  spectral 
peaks  and  possible  vortices  over  the  wing  was  evident,  thus  no  data  ar.e  shown 
in  Figure  152  for  Mach  1.  The  trends  evident  from  Mach  0.7  through  Mach  0.9 
are  thought  to  continue  for  some  increment  in  Mach  number  beyond  0.9,  but  the 
size  of  this  increment  is  unknown.  This  lack  of  data  at  Mach  1.0  is  not 
thought  to  be  critical  because  low  values  of  A^prms  are  indicated  in  Figure 
147.  The  plots  versus  configuration,  which  are  essentially  plots  versus  flap 
deflection,  show  that  at  Mach  0.7  and  in  the  neighborhood  of  12°  wing  angle 
of  attack  the  frequencies  are  relatively  independent  of  flap  position.  The 
snag  vortex  frequency  is  raised  about  fifty  percent  when  the  leading  edge  flaps 
are  lowered  8°.  Values  of  kQ  for  particular  flight  conditions  and  flap 
deflections  can  be  obtained  by  interpolating  among  graphs  in  Figure  152. 

A  wing  other  than  that  of  the  F-4  may  or  may  not  generate  all  three 
vortices  discussed  above.  Most  modern  fighter  wings  will  have  a  tip  vortex 
and  a  leading  edge  vortex.  They  will  probably  not  have  what  has  been  referred 
to  above  as  a  snag  vortex  unless  they  do  in  fact  have  a  snag  or  similar  lead¬ 
ing  edge  discontinuity.  Figures  62,  65,  70,  74,  82,  87,  92  and  96,  and  Figure 
152  indicate  that  the  frequencies  associated  with  the  leading  edge  and  tip  vor¬ 
tices  on  the  F-4  wing  are  essentially  the  same.  It  is  assumed  that  this  re¬ 
lationship  would  hold  true  on  other  fighter  wings.  Thus,  if  a  wing  has  a  lead¬ 
ing  edge  vortex  and  a  tip  vortex  but  no  snag  vortex,  then  for  purposes  of 
modeling,  the  effects  of  the  tip  vortex  can  be  accounted  for  by  the  leading 
edge  vortex;  therefore,  only  the  leading  edge  vortex  need  be  considered. 

5.2.4  Preliminary  Design  Data  -  In  the  preliminary  design  stage  of  a 
fighter  aircraft,  there  are  at  least  two  ways  to  estimate  the  areas  of  the 
wing  over- whiqji  the  fluctuating  pressures  are  affected  by  vortices  and  which 
types  of  vortices  they  are.  The  quickest  way  is  to  use  engineering  Judge¬ 
ment  and  consider  the  characteristics  of  the  wing  (e.g.,  does  it  have  a  snag) 
and  the  regions  of  the  F-4  affected  by  vortex  frequencies  shown  in  Figures 
62,  65,  70,  74,  82,  87,  92  and  96.  For  example,  on  a  wing  with  a  planform 
similar  to  that  of  the  F-4  wing,  but  having  a  smooth  loading  edge,  the  Judge¬ 
ment  for  Mach  0.7  and  aw  -  12°  might  be  as  shown  in  the  following  sketch: 


4l 


Fluctuating  Pressures 
Not  Affected  by  Leading 
Edge  Vortex 


The  second  way  to  estimate  the  types  of  vortices  present  over  a  wing 
and  the  areas  of  the  wing  over  which  they  affect  the  fluctuating  pressures 

is  to  subject  a  model  of  the  wing  to  a  flow  visualization  test.  As  was 

indicated  at  the  beginning  of  this  section,  a  flow  visualization  test  gives 
an  indication  of  the  direction  and  turbulence  of  the  flow.  Interpreting 

this  in  terras  of  the  causes,  whether  vortices  or  other  phenomena,  is  often 

difficult.  It  may  or  may  not  yield  better  results  than  engineering  Judge¬ 
ment  alone, 

5.2.5  Experimental  Values  of  n  -  From  the  plots  of  n  versus  percent 

chord  in  Figures  97  through  1C&,  a  simplified  model  for  the  variation  of  n 

over  the  wing  was  developed  and  is  shown  in  Figure  153.  This  model  is 

similar  to  that  for  the  rms  pressure  distribution  over  the  wing  except  that 

the  values  of  n  at  the  corners  of  the  wing  planforra  are  treated  as  variables 

and  plotted  versus  wing  angle  of  attack,  Mach  number,  and  configuration. 

The  plots  for  configuration  1  show  that  n  =  1  at  the  wing  root  for  all  a 

w 

from  8°  through  lk°  and  all  Mach  numbers  from  0.7  through  0.9.  At  the  wing 
tip  at  Mach  0.9,  *1  is  higher  at  ay  «  12°  than  at  either  8°  or  lh°.  At  the 
tip  at  12°  wing  angle  of  attack,  n  increases  at  the  leading  edge  and 
decreases  at  the  trailing  edge  as  the  Mach  number  increases  from  0.7  through 
0.9.  As  discussed  in  Bection  5.2.3  in  connection  with  kQ,  data  on  spectral 
peaks  at  M  =  1.0  are  not  presented  because  no  pattern  was  evident.  The  trends 
shown  for  n  in  Figure  153  from  Mach  0.7  to  Mach  0.9  can  be  expected  to  continue 
for  some  increment  beyond  M  =  0.9i  but,  again,  the  size  of  the  increment  is  un- 


known.  The  plot?  for  configuration  5  in  Figure  153 ("b )  show  that  at  the  wing 

root  and  the  tip  leading  edge,  n  =  1  for  all  a  from  11°  through  lk°  and  all 

Mach  numbers  from  0.7  through  0.9.  At  the  tip  trailing  edge  at  Mach  0.9,  n 

is  higher  at  =  11°  than  at  l4°.  At  the  tip  trailing  edge  at  a •,  =  11°,  n 

decreases  as  the  Mach  number  increases  from  0.7  through  0.9.  Figure  153(c) 

shows  the  effect  of  model  configuration  on  n  at  Mach  0.7  and  a  about  12°. 

w 

At  the  root  and  leading  edge  of  the  tip,  n  is  the  same  for  configurations  1, 
5,  and  6.  At  the  trailing  edge  of  the  tip  n  is  the  same  for  configurations 
1  and  5  but  lower  for  configuration  6.  Values  of  n  for  particular  flight 
conditions  and  flap  deflections  can  be  obtained  by  interpolating  among 
graphs  in  Figure  153. 

5.3  Spatial  Correlation  Effects 

Conceptually,  the  cross-correlation  between  the  buffet  pressures  at  two 
separated  points  may  be  envisioned  to  embody  a  spatial  decay  effect  and 
possibly  some  convection  effect.  On  this  basis,  a  model  for  the  cross¬ 
correlation  function  may  be  written  in  the  form 


K12  (r,x)  =  a 


°2e 


-^l  e-8^  p(i  -  x ’ ) 


(5-7) 


where  and  denote  the  respective  rms  values  at  two  points ,  Ax  and  Ay 
represent  the  x  and  y  separation  distances  between  the  points,  t '  relates  to 
the  tin©  of  convection  of  the  pressures  from  one  point  to  the  other,  and  p 
denotes  the  normalized  autocorrelation  function.  The  following  is  another 
plausible  form  that  combines  the  x  and  y  dependency  into  a  single  spatial 
function  r, 

R12  (r,t)  =  p(t  -  t')  (5-8) 


where  r  is  simply  the  distance  between  the  two  points. 
The  cross  spectrum  indicated  by  equation  (5-8)  is 
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vhi ch  indicates  a  magnitude 


S12  (r,f) |  =  o1  o. 
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(5-10) 


and  a  phase  angle  0  between  the  Quid  and  Co  power  which  is 


0  =  2TTfT* 


(5-11) 


In  Equations  (5-9)  and  (5-10)  $  is  the  normalized  power  spectral  density- 
function.  Note  that  if  x*  is  a  constant  then  the  phase  versus  frequency 
curve  is  a  straight  line  passing  through  the  origin  with  slope  equal  to 
2irx'.  Note  also  that 

p12  (r.T.,  i  £!k^>  .  e-{IH  =  4  y12(r,f)  (W*| 

12  12 
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where  p  (r,x)  is  the  cross  correlation  of  the  pressures  at  points  1  and  2 
normalized  by  dividing  by  the  product  of  the  ms  values  of  the  pressures 
and  y  (r,f)  is  the  magnitude  of  the  cross  power  spectral  density  normalized 
by  dividing  by  the  square  root  of  +he  product  of  the  power  spectral 
densities  of  the  pressures. 

The  experimental  results  tended  to  fit  the  form  indicated  by 
equation  (5-8).  Attempts  were  therefore  made  to  establish  apparent  convec¬ 
tive  effects  by  considering  time  delays  t '  between  various  transducers. 
Theoretically,  if  an  ordered  convection  field  directed  along  a  line  between 
two  transducers  were  involved,  t*  would  lead  to  a  consistent  convection 
velocity  given  by  Ve  =  r/t ' ,  where  r  is  the  distance  between  the  two  points 
under  consideration.  Studies  using  the  r  and  t*  associated  with  many  pairs 
of  transducers  proved  to  be  unproductive  in  defining  convection  velocities. 

A  more  refined  study,  using  data  from  sets  of  three  transducers,  gave  prom¬ 
ise  Of  defining  local  convection  velocity.  However,  this  study'  was  not 
completed  due  to  the  inconsistent  nature  of  the  flow  field.  This  latter 
study  assumed  that  disturbances  propagate  across  the  transducer  set  in  a 
line,  or  front.  The  results  indicated  that  the  idea  of  an  orderly  consis¬ 
tent  convection  did  not  fit  the  data  for  the  F-h  model  wing.  Apparently  the 
various  vortices,  that  emanate  from  the  leading  edge  mix .In  such  a  hetero- 
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geneous  way  to  make  the  flow  appear  chaotic,.  Thus  the  best  conclusion  that 
can  be  drawn  from  the  test  results  is  that  detailed  consideration  of  con¬ 
vection  cannot  be  made,  and  it  is  therefore  neglected  in  this  study.  Essen¬ 
tially  this  amounts  to  setting  t*  =  0  in  equations  beginning  with  (5-8) • 

Tests  on  other  wings  are  needed  to  determine  if  similar  results  are  found. 

Results  for  spatial  decay  effects  were  found  to  be  somewhat  better, 
although  still  quite  rough.  In  equation  (5-12)  the  spatial  decay  factor 
represents  essentially  a  cross-correlation  coefficient  between  the  two 
pressures,  or 


12  t  ,x  —12  /  s  -<$|r| 

)  (r,x’)  =  y  (r,f  )  =  e  11 


(5-13) 


where  6  is  a  spatial  decay  constant  to  be  empirically  established.  In  this 
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study,  the  average  value  y  of  the  magnitude  of  the  normalized  cross 
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power  spectral  density  was  used.  Values  of  y  derived  from  several 
normalized  cross  power  spectral  density  functions  are  shown  in  Figure  151* 
in  terms  of  the  separation  distance  r.  The  curve  shown  in  this  figure  is 
on  exponential  fit  of  the  data  given  by 


Y12  .  e-°*25r 
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(Although  there  is  scatter  in  the  plotted  data  of  Figure  15  ^ ,  an  exponential 
curve  fit  is  the  most  tractable  representation  of  the  dnta,and  it  has  been 
used  successfully  by  previous  investigators.)  This -expression  Gan  be 
converted  to  one  of  general  use  by  expressing  it  in  non dimensional  fora 
through  use  of  the  mean  aerodynamic  chord,  thus 
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Fbr  the  model  e  «  19.25  inches,  and  thus  in  the  fora  for  general  use  y 
is  given  by 
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Within  the  uncertainty  in  the  value  of  the  decay  factor  6  =  as 

indicated  by  the  scatter  in  the  plotted  points  in  Figure  15 1*,  a  and  6 

h  8 

(see  equation  (5-7)  for  example)  may  also  each  be  set  equal  to  — .  Then 
an  alternative  expression  is 


-12  =  -U.C|  Ax|/c  -4.8|  Ay  |  /  e 
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(5-16) 


Equations  (5-15)  and  (5-16)  represent  different  forms  of  the  spatial  decay 
factor  that  has  been  derived  for  general  use  from  the.  F-l*  tests.  It  is 
desirable  to  determine  if  tests  on  other  wings  yield  a  coefficient  similar 
to  the  U.8  value  derived.  In  application  to  other  wings,  c  would  be  the 
value  of  the  mean  aerodynamic  chord  for  the  wing  under  consideration. 

The  above  results  for  the  cross  correlation  and  cross  power  spectral 
density  functions  for  the  pressures  at  points  1  and  2  separated  by  a 
distance  r  may  now  be  expressed  in  the  terms  of  equations  (5-3)  and  (5-5) 
as  follows 
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where  a.  and  o„  refer  to  the  respective  nns  values  of  pressure  at  the  two 
points,  fi  =  and  the  'ther  parameters  are  as  defined  in  Section  5.2.2. 
5 . b  Procedure  for  Applying  Prediction  Method 

Following  is  a  step  by  step  procedure  to  be  followed  to  predict  the 
fluctuating  pressures  on  a  fighter  airplane  wing  at  a  given  point  in  the 
transonic  speed,  high  angle  of  attack  flight  regime.  The  procedure  is 
based  on  the  prediction  method  described  in  this  section. 


5.^.1  Spectral  81  wipe 
STEP  1 

V  find  c  are  known. 
u'i'ii'i1  y 


Estimate  the  types  and  locutions  of  the  vortices.  To  do  this, 


1*6 


use: 


a)  Engineering  Judgment  considering  the  characteristics  of  the 
wing  (is  there  a  snag?)  and  the  regions  of  the  F-4  wing 
affected  Ly  vertex  frequencies  shown  in  Figures  62,  65,  70, 

74,  82,  87  ,  9‘2.  and  96. 

b)  Flow  visualization  tests  combined  with  engineering  judgement. 

If  only  leading  edge  and  tip  vortices  are  present,  the  tip 

vortex  can  be  neglected.  If  a  snag  vortex  overlaps  another  vortex, 
the  effect  of  the  other  vortex  should  be  neglected  in  the  area  of 
overlap. 

STEP  3 

Establish  the  nondimehsional  vortex  frequency,  k  ,  from  Figure  152. 
Assume  kQ  established  by  a  vortex  to  be  applicable  over  a  chord 
line.  If  a  vortex  does  not  exist  over  part  of  the  chord  line,  kQ 
will  serve  as  a  cutoff  frequency  in  the  model  for  the 
spectral  shape  to  limit  the  frequency  range  over  which  most  of 
the  energy  is  distributed. 

STEP  4 

Establish  n,  the  ratio  of  the  peak  to  the  low  frequency  level 
of  the  spectrum,  from  Figure  153. 

STEP  5 

With  n,  obtain  the  parameter  X  from  Figure  ,150. 

STEP  6 

With  X,  determine  the  spectral  shape  of  the  kQ~-j”  function  from 

Figure  149. 


STEP  7 


With  kQ,  compute  the  nendimensiona't  power  spectral  density  as  n 
function  of  the  nondimensional  frequency  k  =>  — -  : 

k  a  k  (p~) 

0  i 

V 

SUlL) 
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With  e  mid  V,  calculate  the  nomnlined  power  spectral  density 

-44-i-  as  a  function  of  the  frequency  f. 
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SECTION  6 

PREDICTION  METHOD  FOR  BUFFET  DYNAMIC  RESPONSE 


The  structural  dynamic  response  prediction  technique  chosen  for  this 
first  generation  buffet  prediction  method  is  based  on  work  that  is  reported 
in  References  10  and  21.  That  analysis,  derived  to  estimate  the  acceleration 
response  to  buffet  pressures  of  truncated  conical  shells,  incorporates  approx¬ 
imations  that  are  considered  consistent  with  the  goals  of  this  development. 

The  problem  formulation,  mathematical  modeling,  and  simplifying  assumptions 
are  discussed  in  this  section,  and  equations  of  motion  are  presented.  It  is 
to  be  noted  that  essentially  all  of  the  development  of  the  mathematical 
response  equations  presented  in  this  section  was  done  in  the  above  mentioned 
References  10  and  21.  Additional  information  is  given  in  Reference  22. 

6.1  Outline  of  Response  Prediction  Method 

The  steps  involved  in  the  response  analysis  are: 

(1)  Divide  the  wing  planform  into  a  gridwork  of  surface  areas  (see 
Paragraph  6,2). 

(2)  Define  structural  properties  of  the  aircraft.  The  elastic  vibration 
characteristics  of  the  aircraft  wing  are  specified  in  terms  of 
uncoupled  modes.  The  required  characteristics  are  modeshapes 

^  (x,y),  modal  frequencies  generalised  masses  f*L,  and  modal 
damping  coefficients  1^.  These  data,  most  conveniently  defined 
for  the  case  of  zero  flight  speed,  may  be  obtained  either  experi¬ 
mentally  (for  an  existing  airplane)  or  analytically. 

(3)  Define  the  fluctuating  pressures  within  each  area  element  in  terms 

of  the  pressure  power  spectra  (specified  as  normalized  power  spectral 
densities,  overall  rmo  fluctuating  pressure  coefficients,  freestream 
dynamic  pressure,  and  spatial  decay  constant). 

(I*)  Determine  generalised  aerodynamic  forces  using  pressure,  area  and 
modes hape  definitions.  The  essentials  of  the  mathematical  treatment 
used  in  the  response  prediction  method  are  presented  below  in 
Paragraph  6.3. 

(5)  Compute  response  power  spectra  and  ms  values  of  response 

parameters  at  selected  locations  on  the  structure.  Expressions  for 
the  response  spectra  are  discussed  in  Paragraph  6,3, 
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6.2  Wing  Area  Sectioning 

The  wing  planform  is  divided  into  a  number  of  basic  areas  for  the  pur¬ 
pose  of  simplifying  the  computation  of  generalized  aerodynamic  forces.  While 
the  area  sectioning  may  be  done  in  several  ways,,  the  method  selected  for 
this  study  defines  the  boundaries  of  the  areas  with  a  series  of  constant  per¬ 
cent'  chord  lines  intersecting  with  streamwise  lines  (See  Figure  155).  The 
basic  areas,  indicated  by  numbers  at  their  centers,  are  also  members:  of  larger 
surface  areas  (designated  as  regions )that  are  indicated  by  numbers  in  circles. 
The  regions  are  areas  over  which  the  spatial  decay  coefficients  and  the  shapes 
of  the  power  spectral  densities  of  pressure,  both  discussed  in  Section  5,  are 
assumed  to  be  constant.  The  basic  areas  are  sized  so  that  cross  correlation 

— 8  i*  /(j 

between  these  areas  can  be  neglected.  A  decay  value  of  e  *  <_0.3  gives 

a  value  of  r  approximately  0.25c  so  that  basic  areas  with  dimensions  of  about 
one-fourth  the  mean  aerodynamic  chord  should  fulfill  this  latter  requirement. 
It  should  be  noted  that  gross  structural  response,  and  not  local  response, 
is  of  most  concern  so  that  localized  effects  (such  as  correlation  immediately 
across  a  basic  area  boundary)  will  not  significantly  affect  overall  wing 
response. 

6.3  Overview  of  Mathematical  Development 

The  aircraft  structure  has  normal  modes  ^(Xjy)  corresponding  to  a  set 
of  normal  coordinates  q^(t),  The  modal  response  due  to  pressure  excitation 
p(x,y,t)  is  determined  from  the  fundamental  equation  in  the  time  domain, 


q1(t)  +  2Z1  wi  4t(t)  +  wi2  qi(t)  =  If  4^(x,y)  p(x,y,t)  dxdy  (6-1) 

wing 

planform 

where 

u»i  =  frequency  of  mode 

M,,  ®  generalised  mass  of  mode 

A  1 

Cj  a  damping  ratio  of  mode  4^ 

Quae  the  modal  responses  q^(t )  are  determined,  then  the  local  displacement 
response  s  (t)  at  any  point  'V  on  the  wing  is  determined  by 


s(4(t)  al^ai) 


(6-2) 
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which  states  that  the  total  response  at  any  point  "a"  is  the  sum  of 

the  products  of  modeshape  values  at  "a”  times  the  corresponding  generalized 

coordinate  values  obtained  from  (6-1). 

Although  equations  (6-1)  and  (6-2)  are  presented  in  the  time  domain,  the 
random  nature  of  the  variables  makes  the  solution  of  these  equations  not 
adaptable  to  deterministic  techniques.  Rather,  the  nature  of  the  measured 
buffet  pressures  indicates  that  statistical  methods  of  analysis  are  required. 
The  well  known  procedures  of  power  spectrum  analysis  have  been  used  in 
References  10  and  21,  and  are  followed  here.  The  resulting  equations, 
expressed  in  the  frequency  domain,  are  obtained  by  applying  conventional 
Fourier  transform  techniques  to  the  time  functions.  Therefore,  the  relation¬ 
ships  that  follow  are  presented,  in  the  form  of  power  pectrum  functions  that 
are  derived  from  the  conventional  time  varying  functions. 

The  following  relationship  may  be  considered  as  a  frequency  domain  counter¬ 
part  of  (6-1)  because  it  relates  modal  spectrum  response  to  pressure  spectrum 
forcing  functions: 


S*J(w)  = 
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Ilf  f  IS  (x,y,x  ,y  ,u>)  <i>.  (x,y) 


i  i 


•  <Mx  ,y  )  dA^  dA  % 
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where 


i  1 

S^J(w)  =  the  response  cross  spectrum  for  the  normal  coordinates  q^  and 

A 

gtt(x,y,x  *y  i^)  =  pressure  cross  spectrum  between'  points  x,y  and 

x-.y- 

and  the  %  are  frequency  response  functions  for  modes  i  and  J , 
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..‘tug  areas  are  designated  by  A^.  and  A,,  and  the  indices  k  and  &  both  vary  from 
1  to  the  total  number  of  basic  areas .  The  frequency  domain  counterpart 
-f  equation  (6-2)  is  expressed  as 


51 


••  -'S'  ' ‘■'f'  •■*’■ ''  '•  ■■ r’  -:  '■'■,  •  ■  :  ■;•■  ■  ■  ’  ■ . ,  ■■’ ; '■ 

****<0*— ..-«.)  -  -  -  ■-  --  -  -  -  ■  ■  ■  -■  —  ■•■■■•■■■'  ■  . - . -  .  ~ _ „.  . . •„•, r,.nv^w,..>.i,w.i. 


ab 


i«5 


Sz  U)  =  U^j  .[sq 


(«)]. 


lV 


T 


(6-6) 


in  which 


4>  .  =  the  ith  degree  of  freedom  modeshape  value  at  point  "a” 
ax 


,ab/ 


and  S  (to)  =  displacement  response  cross  spectrum  between  points  a  and  b. 
z 


6.3.1  Modeshape  Averaging  -  A  simplication  to  the  form  of  (6-3)  is 
introduced  by  a  modeshape  averaging  technique  that  is  introduced  in 
Reference  10.  In  this  method,  it  is  simply  assumed  that  a  modeshape  dis¬ 
placement  value  is  constant  within  any  basic  area  A^.  The  value  of  the 
modeshape  used  for  any  mode  i  within  area  A^  is  taken  as  an  average  value 
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Since  the  modeshape  values  within  an  area  then  become  independent  of  x  and  y, 
equation  (6-3)  can  be  rewritten 
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With  the  following  definition 
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equation  (6-7)  can  be  expressed  in  the  fom 
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6.3.?.  Real  Representation  of  Response  Spectrum  -  Reviewing  equations 
(6-6)  and  (6-9), which  give  the  deflection  response  cross  spectrum  between 
pointo  a  and  b,  it  should  be  recalled  that  the  cross  spectrum  is  in  general 
a  complex  quantity.  Therefore,  the  cross  spectrum  of  any  quantity,  whether 
>  insure  or  displacement,  can  be  written  in  terms  of  a  real  and  an  imaginary 
part,  vis,. 


:*U(«)  *  C**iu)  -  i^(oi) 
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The  main  interest  here  is  the  power  spectrum  of  the  response  at  a  point, 
which  is  obtained  by  letting,  for  example,  b  =  a  in  the  above  equations.. 

By  definition  a  power  spectrum  is  a  real  function;  hence  only  the  real  part  . 
of  equation  (6-6)  is  considered  in  what  follows.  Since  the  only  imaginary 
terms  in  equation  (6-6)  are  contained  within  [SX^(u)],  the  real  part  of  this 
equation  for  the  special  case  of  b  =  a  can  be  written  as 

c^U)  =  U^J  [cTJU)  L <j>aJ JT  "  ’(6—11) 

Measured  spectra  are  expressed  only  for  positive  frequencies,  and  are 
usually  presented  in  terms  of  frequency  in  Hertz.  Therefore  it  is  expedient 
at  this  point  to  define  the  spectrums  for  positive  frequencies,  and  for 
frequency  in  Hertz.  Thus  equation  (6-11)  becomes 

c“(f)  =  U^J  [c*J(f)3  (6-12) 


Following  the  development  in  Reference  21,  the  real  part  of  equation 
(6-9)  can  then  be  expressed  in  real  terms,  viz., 
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+  Hg(f)  Z  l  Tik  ^  S  f  (x,y,x'y',f)  dA^  dA£]  (6-13) 
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where 
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(Note:  Hg( f)  =  0  when  i  «  J.) 
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Equation  (6-13)  is  an  expanded  form  of  equation  (6-3)  making  it  more  adaptable 


to  the  use  of  measured  data 
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6.3.3  Simplifying  Assumption  — It  is  now  assumed  that  the  cross 
correlation  function  (or  normalized  cross  power  spectral  density  magnitude) 
of  the  fluctuating  pressures  is  negligibly  -small  except  when  l  =  k.  This 
is  the  case  if  the  areas  are  chosen  sufficiently  large.  The  rationale  used 
for  the  selection  of  area  size  based  on  decay  data  was  presented  in  Para¬ 
graph  6.2.  The  cross  correlation  function  between  any  two  points  lying 
within  a  basic  area  is  represented  as  the  product  of  the  autocorrelation 
function  at  the  center  of  A^  and  exponential  decay  functions  of  the  fom 
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The  cross-correlation  can  then  be  written  as 


Rp  (y,v,r)  =  Rq(t)  e"W  e“^ 


(6-17) 


where 


y  =  j  x— x’ 
\j  =  ly— y  * 


(6-18) 
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and  Rq( t )  is  the  autocorrelation  function  at  the  center  of  A^. 

The  equation  for  the  cross  power  spectrum  given  in  equation  (3-6)  of 


Section  3.3  can  be  written  as 


S  (u,v,f)  a  /  R  (y,v,r)  e”i27rfl  a* 


o  /  R  (y,v,r)  (co32nft  -  i  sin2nfT)  dr 
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(6-20) 


The  real  and  imaginary  parts  of  this  equation  are  the  Cp(x,y,x' ,y*  ,f)  and 
the  Qp(x,y*x' ,y' ,f) ,  respectively ,  of  equation  (6-13). 

Substituting  equation  (6—17)  into  equation  (6-20),  integrating,  and 
setting  £  =  k  gives 
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and  Qp(x,y,x  ,y  ,  f)  =  0 


(6-22) 


where  Cp  (f)  is  the  power  spectrum  at  the  center  of  A^,.  The  function 
Qp  (x,y,x  ,y  ,  f)  is  zero  because  it  comes  from  the  integration  of  the 
product  of  an  even  function,  Rq(t),  times  an  odd  function,  sin  27rfT,  over  the 
range  (-»,“). 

Using  equation  (6-21)  in  the  first  integral  of  equation  (6-13)  yields 
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The  integral  on  the  right  side  of  this  equation  is  called  the  decay 
correction  factor  for  A^,  DCFk> 
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To  simplify  the  integration,  A^  is  assumed  to  be  rectangular  with  chordwise 
dimension  Ax  and  spanwise  dimension  Ay.  Ay  is  the  distance  between  the 
constant  percent  span  lines  which  bound  A^,  Ax  is  an  average  distance 
between  the  constant  percent  chord  lines  which  bound  A^. 

The  resulting  integration  yields 
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Equation  (6-13)  can  now  be  written  in  the  form 
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It  is  indicated  in  References  10  and  21  that  a  further  simplification 

is  made  by  neglecting  the  values  of  C^(f)  when  i  4  «J.  Then  equations  (6-12) 

Q. 

and  (6-26)  can  be  written 
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The  computer  program  modified  for  canputing  the  response  spectrums  provides 

i  1 

the  option  of  using  or  not  using  C  J(f)  when  i  4 

<1 

The  acceleration  response  power  spectra  are  found  from  equation  (6-12)  to 
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where 


ciJ(f)  =  (2tt)U  f*  c*J(f) 
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With  the  same  simplifications  used  to  obtain  equation  (6-28),  equation  (6-30) 


can  be  written  as 
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where 
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The  mean  square  responses  are  given  by  the  integrals  ei'  equations  (6-28), 
(6-29),  (6-32),  and  (6-33) 
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SECTION  7 

PREDICTION  OF  AIRCRAFT  RESPONSE  AND  COMPARISON  WITH 
MEASURED  FLIGHT  DATA 


The  method  presented  in  Section  6  was  used  to  predict  the  response  of 
the  F-h  aircraft  to  buffeting  pressures  at  three  flight  conditions.  The 
method  presented  in  Section  5  was  used  to  predict  the  fluctuating  pressures 
which  were  used  as  input  data  in  the  response  predictions.  The  predicted 
responses  were  compared  with  measured  flight  acceleration  data  to  provide 
an  indication  of  the  validity  of  the  buffeting  pressure  and  response  pre¬ 
diction  methods  and  the  results  are  presented  in  this  Section. 

7.1  Response  Equations 

Equations  (6-28),  (6-29),  (6-32),  and  (6-33)  were  used  to  compute  the 
response  of  the  F-4  aircraft. 

7.2  Wing  Idealization 

The  F-ii  wing  was  idealized  as  having  a  straight  leading  edge  and  being 
divided  into  3  regions  and  l8  areas  as  shown  in  Figure  156.  The  power 
spectral  density  shapes  in  regions  1  and  2  are  determined  by  the  leading 
edge  and  snog  vortices,  respectively ;  and  the  shape  in  region  3  is  deter¬ 
mined  by  either  the  tip  or  snag  vortex,  depending  on  their  relative 
dominance.  The  gridwork  of  areas  is  a  compromise  that  provides  reasonable 
definition  of  the  structural  vibration  mode  shapes  and  the  overall  rms 
fluctuating  pressure  distribution  while  limiting  the  maximum  value  of  the 
normalized  cross-correlation  function  for  tire  pressures  at  the  center's  of 
adjacent  areas  to  about  0.2  inboard  and  0.6  outboard. 

7.3  Structural  Representation  of  F-i>  Wing 

Data  for'  the  first  10  symmetric  and  the  first  10  antisymmetric  normal 
modes  of  the  F-l*  wing  with  full  wing  fuel  were  used  to  represent  the 
structure  responding  to  the  random  excitation*  The  data  are  based  on  a 
finite  element  analysis  program  being  developed  for  the  wing  in  con¬ 
junction  with  the  production  F-i*  effort,  mid  they  represent  the  best 
analytical  model  of  the  wing  to  date.  The  first  four  analytical  and  mea¬ 
sured  frequencies  are  shown  in  Table  X  and  are  in  good  agreement. 

Table  XT  shows  the  natural  frequencies  and  generalised  masses  for  all 
m  modes  ordered  mi  used  in  this  study.  Table  XII  present*  the  mode- 
shapes  for  these  moden  in  a  column  format.  The  column  numbers  in  Table 
XII  correspond  to  the  mode  numbers  in  Table  XI,  Table  XIII  relates  the 


#**  ' 


Table  X 

Comparison  of  Analytical  and  Measured  F4  Wing  Vibration  Frequencies.  H  z 


Wing 

Mode 

Predominant 

Characteristic 

Full  Wing  Fuel 

Symmetric 

Antisymmetric  | 

F4H-1 

Test 

F-4E 

Test 

F-4E 

Theory 

F4H-1 

Test 

F-4E 

Test 

F-4E 

Theory 

1 

IstBending 

8.05 

8.3 

8.4 

10.5 

10.0 

* 

10.4 

2nd  Bending 

19.2 

17.1 

17.5 

20.3 

17.8 

19.9 

1st  Torsion 

25.6 

28.8 

29.9 

28.2 

30.0 

30.5 

Hi 

3rd  Bending 

- 

35.0 

35.7 

- 

- 

38.3 

TaoioZI 


Analytical  F-4  Wing  Vibration  Frequencies 
and  Corresponding  Generalized  Masses 


Mode 

No. 

Symmetric  or 
Antisymmetric 

— 

Natural 

hequency 

(Hz) 

Generalized 

Mass 

^Lb-Sec2j 

1 

Symmetric 

8.38 

0.299 

2 

17.5 

0.218 

3 

•J9.9 

0.372 

4 

35.7 

0.118 

5 

47.8 

0.101 

6 

57.3 

0.256 

7 

62.6 

0.338 

8 

65.8 

0.0888 

g 

79.6 

0.186 

10 

94.0 

0.168 

ii 

Antisymmetric 

10.4 

0.218 

12 

19,9 

0.319  - 

13 

30.5 

0.202 

14 

38.3 

0.134 

15 

48.7 

0.110 

16 

57.8 

0.282 

17 

64.1 

0.343 

18 

66,0 

0.0991 

19 

79.8 

0.189 

20 

94.7 

0.177 
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Null's. 

1.  Voi  ticai  1  lanslationx  .in'  I'nsilivti  I1f>  (Rajhi  Winn),  Pitch  Rotation  is  Positive  Nosu 
Up,  anil  Holl  Hmnlion  is  Potihvo  R19IH  Wnnj  lip  Up. 

'2.  Values  Shown  (7)1  Veriiral  Tiunsl.ninns  of  Points  on  tha  Wino  I  minute  lit  touts  ot 
LsoPi  Wn>ii  Ewsiir  I'.Vtoi millions  an  J  Rim<!  Airpl.mn  Motions, 
rim  uilnnms  liuhr-tto  Mm  In  Nttmlwi  anil  thu  Rows  liuheatu  Alrplano  Coordinate. 
4,  P.i'o  Tallin  H  1 1  Ini  ilia  Nut ural  t-rui|u,-noy  anil  Gunoruli/od  Mass  of  (inch  Moil*. 

I'm*  I'aPIn  H  IV  tor  How  litimtitioulion. 

I!  Kxiioiiu.iiiol  Nouintin  is  nsod,  ./.4470-03  ”  -7.447  x  10  ^ 
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Table  Xm 

Identification  of  Rows  in 
Modeshape  Matrix  (Table  XU) 


Row  of  Matrix  Airplane  Coordinate 

1  Vertical  T ranslation  of  Center  of  Area  1 

2  2 

3  3 

4  4 

5  5 

6  6 

7  7 

8  8 

9  9 

10  10 

11  11 

12  12 

13  13 

14  14 

15  15 

16  16 

17  17 

18  18 

19  Vertical  Translation  of  Flight  Accelerometer 

20  Rigid  Airplane  Vertical  Translation 

21  "  "  Pitch  Rotation 

22  '*  “  Roll  Rotation 


rovs  of  Table  MI  (i.e.,  the  rows  of  the  modeshape  matrix)  to  airplane 
coordinates.  The  definition  of  the  moaeshapes  at  these  coordinates  (l8 
area  centers  and  the  flight  accelerometer)  was  obtained  by  linear  inter¬ 
polation  among  the  original  57  wing  coordinates  in  the  vibration  solution, 
and  the  values  shown  were  rounded  to  3  significant  figures  in  the  response 
computations. 

Structural  damping  coefficient  (g^j  data,  from  Reference  23,  are 
available  for  only  the  first  two  or  three  symmetric  and  antisymmetric  wing 
vibration  modes  and  are  presented  in  Table  XIV.  Since  the  measured 
structural  damping  values  increase  with  mode  number  through  the  third 
symmetric  and  second  antisymmetric  modes  and  since  they  cannot  increase 
without  limit,  structural  damping  coefficients  for  the  higher  modes  have 
been  estimated  as  0.1,  which  is  approximately  the  largest  value  measured 
for  the  first  three-  mode's-; 

The  total  damping  for  an  airplane  in  flight  is  due  to  both  structural 
damping  and  aerodynamic  damping.  Aerodynamic  damping  coefficients  (g  ) 

ft 

for  the  F— U  wing  vibration  modes  used  in  this  study  were  obtained  from 
velocity  versus  damping  flutter  solutions,  and  values  for  the  three 
dynamic  pressures  investigated  herein  are  given  in  Table  XIV.  Also  shown 
in  Table  XIV  are  the  total  modal  dumping  ratios 


Structural  damping  in  larger  than  aerodynamic  damping  in  all  but  the  first 
one  or  two  symmetric  and  antisymmetric  modes.  Response  calculations  were 
also  made  using  the  name  value  of  c,  in  all  modes  because  there  is  consider¬ 
able  uncertainty  in  the  damping  values  shown  In  Table  XIV  and  there  would 
be  much  more  uncertainty  In  the  preliminary  design  stage  of  a  new  airplane, 
v.  It  1  i  i  ghl  don 'H  1  t '  ns 

The  response  of  the  i-'J»  wing  to  buffeting  pressures  was  predicted  for 
<  ho  I'hrr  -  .  Light  f  it  or.-;  .".hewn  y.i  ‘•‘■ih.lc  XV,  These  n  ight  conditions  are 
ones  for  which  measured  acceleration  power  spectral  density  data  were 
obtained  n'  the  flight  accelerometer  in  the  flight  tost  program  reported 
in  Reference  2.  For  purposes  of  predicting  the  response,  the  first  two 
flight  conditions  were  considered  to  have  the  same  Mach  number  (0.8)  and 
While  of  attack  (I0°i.  They  are  conditions  of  moderate  buffet.  The  third 
is  a  condition  of  heavy  buffet. 


Tab 

F4  Wind  Moda 


Frequency  (Hzl 

Structural  Damping 
Coefficient,  gs 
(Full  Wing  Fuel)  <1> 

Measured 

jAnalytical 

8.05 

8.38 

0.0172 

19.2 

17.6 

0.0- 

1 

25.6 

29.9 

0.084 

35.7 

0.1 

(2) 

47.8 

67.3 

62.8 

66.8 

79.6 

94.0 

10.5 

10.4 

0.030 

20.2 

19.9 

0.108 

30.5 

0.1 

(2) 

38.3 

48.7 

57.8 

64.1 

66.0 

79.8 

94.7 

Aerodynamic  Damping 
Coefficient,  ga 


9s  +  9a 


q  =  470 
psf 


q  =  220 
psf 


S! 


7. 5  Prei 

Tlie  prrcedures  of  Section  5,4  were  used  to  predict  the  normalized 

G( 

power  spectral  density  ■  - for  each  region  and  ACL.  for  each  area. 

a2  rrms 

The  boundaries  of  the  leading  edge,  tip,  and  snag  vortices  were  estimated 
based  on  Figures  6 2,  65,  70,  and  7*+.  The  boundaries  for  flight  conditions 
1  and  2  are  shown  in  Figure  157,  and  those  for  flight  condition  3  are  shown 
in  Figure  158.  The  vortex  types  and  corresponding  values  of  kQ,  n,  and  X 
associated  with  each  region  are  shown  in  Table  XVI.  The  values  of  n  and  X 
associated  with  a  region  are  those  evaluated  at  the  center  of  the  region 
since  for  the  structural  res;  use  computations  the  spectral  shape  of  the 
fluctuating  pressures  is  assumed  constant  for  a  region-.  The  values  of 
^^rms  oen'ter  each  area  are  listed  in  Table  XVII.  The  values  are 

t.’ie  same  for  all  three  flight  conditions.  The  mean  aerodynamic  chord  of  the 
F-4  wing  is  192.5  inches.  Hence,  the  streamwise  and  spanwise  decay 
coefficients  are 


„  _  ,  liS.125  „  .  -1 

“  = "  ■ 5  “  iw  “  °-25  ln 


7.6  Results 

7.6.1  Effect  of  limping  -  The  effect  of  the  choice  of  damping  ratio 

C  on  the  response  predictions  is  shown  in  Figure  159.  All  three  curves  are 
predictions  of  acceleration  power  spectral,  density  at  the  flight  aceelero- 
metei  for  flight  condition  1.  In  each  case,  the  damping  ratio  was  the  same 
for  all  20  elastic  wing  modes.  The  values  used  were  0.02,  Q.05,  and  0.1.  The 
curves  are  all  higher  than  the  corresponding  measured  data.  Of  the  three 
curves  shown,  that  for  f,  =  0.05  best  matches  the  shape  of  the  measured 
spectrum  in  that  the  ratios  of  peaks  to  adjacent  valleys  in  the  spectuma  are 
similar.  This  is  shown  in  Figure  161. 

7.6.2  Variations  in  Measured  Response  Data  -  Measured  acceleration 
power  spectral  density  curves  for  the  flight,  accelerometer  at  flight  condi¬ 
tions  1  and  2  are  shown  in  Figure  160.  These  flight  conditions  arc 
essentially  the  same  except  they  are  at  different  altitudes  and  hence  have 
■.:!  fferont  dynamic  pressures.  An  discussed  in  Section  4.H.1,  the  overall  rms 
level  of  the  fluctuating  pressures  is  generally  assumed  to  be  proportional  to 
the  freest ream  dynamic  pressure,  q.  Other  things  being  equal,  the  rms  level 
of  the  response  should  also  be  proportional  to  q.  One  might  thus  expect  the 
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Flight  Conditions  for  Which  Responses  of  F4  Wing 
to  Buffeting  Pressures  Were  Predicted 


Flight 

Condition 


Aircraft 

Configuration 

Mach  No. 

M 

Wing 

Angle 
of  Attack 
(deg) 

Altitude 

(ft) 

Velocity 

(ft/sec) 

Dynamic 

Pressure 

(psf) 

1 

0.82 

9.5 

25,000 

830 

350 

1 

0.79 

10 

35,000 

770 

220 

j 

1 

0.92 

10.5 

26,500 

930 

470 

Constants  Defining  Spectral  Shapes 


Flight 

Condition 


Region  Vortex  Type 


Leading  Edge 
Snag 
Snag 


Leading  Edge 
Snag 
Snag 


1.0  m 
1.6«> 
1.9 


( 1)  Difference  Between  tills  Value  of  a  and  the  Corresponding  Value  for  Flight  Conditions  1  and  2 
is  due  to  Uncertainty  in  Graphical  Interpolation  from  Figure  160. 


1 


Table  23ZII 

ACprm$  at  Centers  of  Areas  for  Flight  Conditions  1, 2  and  3 


Area 

ACd 

— 

prms 

1 

0.052 

2 

0.056 

3 

0.060 

4 

0.056 

S 

0.069 

6 

0.081 

7 

0.060 

8 

0.079 

9 

0.098 

10 

0.062 

11 

0.086 

12 

0.11 

13 

0.064 

14 

0.093 

15 

0.1? 

16 

0.066 

17 

0.098 

18 

0.13 

U 


vW  curves  in  Figure  l60  to  have  the  same  shape,  with  the  overall  rms  for 
flight  condition  1  being  higher  than  that  for  flight  condition  2  by  the 
ratio  of  the  q’s ,  i.e.,  by  350/220  =  1.6.  In  fact,  the  curves  are  similar 
but  not  identical  in  shape,  and  they  have  the  same  overall  rms  level.  While 
the  reason  for  this  is  not  clear,  part  of  the  difference  between  the  actual 
and  expected  relationship  between  the  two  curves  may  be  due  to  aero- 
elastic  effects.  Part  may  be  due  to  the  fact  that  the  Mach  number  and  angle 
of  attack,  while  nominally  0.8  and  10°,  respectively,  for  both  q's  are 
actually  slightly  different.  Another  possibility  is  that  the  rms  fluctuating 
pressure  coefficient  is  not  independent  of  altitude. 

746,3  Comparison  of  Predicted  and  Measured  Response  Data  -  Comparisons 
of  predicted  and  measured  acceleration  response  at  the  flight  accelerometer 
for  flight  conditions  1,  2,  and  3  are  shown  in  Figures  161,  162,  and  163, 
respectively.  Predicted  spectrums  are  shown  for  two  different  modal  damp¬ 
ing  distributions.  One  distribution  was  obtained,  as  discussed  previously, 
from  the  sum  of  the  structural  and  aerodynamic  damping  in  each  mode  and  is 
tabulated  in  Table  XIV.  The  other  is  a  constant  damping  ratio  of  0.05  used 
for  all  20  modes.  The  value  of  0.05  was  chosen  because  it  is  in  the  range 
of  values  that  have  been  found  to  be  representative  and  the  resulting 
ratios  of  peaks  to  adjacent  valleys  agree  well  with  those  for  the  measured 
data. 

Two  observations  regarding  the  comparisons  shown  in  Figures  l6l,  162, 
and  163  are  (l)  the  shapes  of  the  predicted  response  spectrum:)  are  in  rea¬ 
sonable  agreement  with  those  os’  the  measured  response  spectrums  except 
over  the  middle  frequency  band  while  (2)  the  predicted  responses  are  higher 
than  the  respective  measured  responses. 

The  differences  between  the  predicted  and  measureu  spectral  shapes  in  the 
2  to  8,  8  to  12,  20  to  30,  and  30  to  1*0  Hz  frequency  bands  ore  now  discussed 
in  detail.  In  the  2  to  8  Hz  band,  the  predicted  spectrums  decrease  rapidly  with 
decreasing  frequency  because  no  structural  modes  below  8.1*  Hz  were  used  in 
the  response  study.  The  measured  spectrum  for  flight  condition  2  also 
decreases  rapidly  with  decreasing  frequency  below  8  Hz,  but  those  for  flight 
conditions  1  and  3  decrease  slightly  from  8  Hz  and  then  remain  essentially 
flat  down  to  2  Hz,  Possible  contributory  causes  of  tire  latter  behavior 
arc  the  fuselage  first  torsion  mode  (which  has  significant  wing  and  atabilator 
motion)  at  6.3  Hz, rigid  airplane  short  period  and  roll  subsidence  modes 


between  0.3  and  0.8  Hz,  low  frequency  bias  in  the  data,  and  the  smoothing 
el  feet  of  the  4  Hz  bandwidth  used  in  the  spectral  analysis. 

In  the  8  to  12  Hz  band,  the  predicted  spectrums,  particularly  those 
based  on  0.05  damping  ratio,  have  2  peaks  (due  to  symmetric  and  antisymmetric 
wing  first  bending  modes)  while  the  measured  spectrums  have  only  one.  That 
only  one  peak  exists  in  the  measured  spectrums  can  be  explained  by  two  facts. 
The  first  is  that  the  4  Hz  analysis  bandwidth  used  in  computing  the  measured 
spectrums  is  considerably  broader  than  the  two  predicted  peaks.  The  second 
is  that  the  total  damping  ratio  in  these  two  modes  is  considerably  higher 
than  the  0.05  value  used  in  the  constant  damping  ratio  cases. 

The  measured  spectrums  contain  a  peak  in  the  neighborhood  of  24  to  26 
Hz  that  is  not  contained  in  the  predicted  spectrums.  It  is  not  contained  in 
the  predicted  spectrums  because  none  of  the  analytical  elastic  wing  modes 
used  in  the  response  calculations  has  a  corresponding  frequency  between  20 
and  29  Hz.  However,  ground  vibration  test  data  show  that  the  fuselage 
second  vertical  bending  mode  has  a  frequency  of  22  Hz  and  that  it  contains 
large  stabilator  motion  and  some  motion  at  the  wing  location  where  the  mea¬ 
sured  data  shown  here  were  obtained.  Fluctuating  pressures  in  the  wing 
wake  impinging  on  the  stabilator  could  excite  this  mode.  Allowing  for  some 
shift  in  frequency,  this  could  cause  the  peaks  at  24  to  26  Hz  in  the  measured 
response  spectrums.  Ground  vibration  test  data  also  show  that  the  fuselage 
second  torsion  mode  has  a  frequency  of  2 6  Hz.  This  is  closer  to  the  peaks 
in  the  measured  response  spectrums}  but,  although  this  mode  contains  large 
stabilator  motion,  it  has  only  negligible  wing  motion  at  the  flight  accelero¬ 
meter  wing  location,  where  the  measured  data  shown  here  were  obtained. 

The  predicted  acceleration  response  spectrums  contain  a  broad  hump 
composed  of  two  peaks  between  30  olid  40  Hz.  The  symmetric  and  antisymmetric 
wing  first  torsion  and  wing  third  bending  frequencies  are  in  this  band.  The 
measured  response  spectrums  3hown  here  have  only  a  volley  between  30  and  40 
Hz,  However,  environmental  vibration  data  from  another  flight  test 
(Reference  24)  show  that,  during  wind-up  turns  at  about  5000  feet  and  about 
2,j ,000  feet  and  subsonic  and  transonic  Mach  numbers,  there  are  strong  peaks 
,’t  the  range  of  30  to  40  Ha  at  about  mid-chord  by  the  wing  fold  and  at  the 
tip  close  to  the  .leading  edge.  This  indicates  vibration  modes  in  this 
jVequoncy  range  have  strong  responses  which  would  also  be  measured  at  the 
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wing  flight  accelerometer  unless  this  location  is  close  to  node  lines  for  these 
modes.  If  this  location  were  close  tothesenode  lines  on  the  airplane, 
then  the  accelerometer  would  measure  little  or  no  response  in  this 
frequency  range  and  there  would  be  a  valley  in  the  spectrums.  However, 
ground  vibration  test  data  indicate  this  location  is  not  on  node  lines  for 
these  modes  in  still  air.  It  is  possible  these  mode  shapes  change  at  air¬ 
speed;  but  at  these  relatively  low  dynamic  pressures,  it  is  doubtful 
that  they  change  enough  that  the  node  lines  all  cross  the  measurement 
location.  In  short,  the  measured  flight  data  of  Reference  2k  indicate 
strong  wing  response  between  30  and  Uo  Hz,  but  no  completely  satisfactory 
explanation  has  yet  been  found  for  the  absence  of  such  response  in  the 
data  from  the  Reference  2  flight  test  program  presented  here. 

The  predicted  rms  responses  based  on  structural  and  aerodynamic  damping 
are  1*.6,  3.1,  and  1.9  times  the  measured  responses  for  flight  conditions  1, 

2,  and  3,  respectively.  The  predicted  rras  responses  based  on  s  =  0.05  for 
all  modes  are  roughly  20  percent  higher  than  the  corresponding  ms  responses 
based  on  structural  and  aerodynamic  damping.  Part  of  the  reason  that  the 
predicted  responses  are  higher  than  the  measured  responses  is  that  the  model 
for  rms  pressure  coefficient  shown  in  Figure  1^7  was  deliberately  selected 
to  be  somewhat  conservative  in  that  it  is  more  repres entative  of  the  high 


fluctuating  pressure  coefficient  conditions.  Another  part  of  the  reason 
may  be  that  in  the  tunnel ,  on  which  the  model  in  Figure  1^7  is  based, 

is  higher  than  in  flight.  As  shown  in  Figures  133  through  ito,  this  is  true 
at  some  flight  conditions  mid  some  points  on  the  wing.  It  was  indicated  in 
Section  fe.U.l  that  this  may  be  due  only  to  slight  differences  in  the  location 


of  the  snag  vortex.  There  are  insufficient  flight  data  to  be  certain.  If, 

however,  A0«  over  the  entire  wing  is  higher  in  the  tunnel  than  in  flight, 
5  rms 

this  would  help  explain  the  high  predicted  responses. 


7.6,1*  nummary  of  Comparison  with  Measured  Response  Data  -  In  summary. 


the  predicted  response  of  the  F— U  airplane  in  moderate  and  heavy  buffet  has  been 


compared  to  measured  data.  There  is  reasonable  agreement  between  the  predicted 
aii-i  measured  spectral  shapes  except  for  three  or  four  frequency  bands, 
pi  ■*!*>; n-nees  over  these  frequency  bands  are  thought  to  be  possibly  due  to 
data  reduction  and  to  differences  between  the  aircraft  in  flight  and  its 


structural  representation  in  this  study.  Although  the  predicted  overall  ms 
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level  is  higher  than  that  measured,  it  is  in  quite  reasonable  agreement 
for  the  more  severe  buffet  conditions.  Some  questions  remain  regarding 
the  effects  of  altitude  on  the  rms  response  and  the  effects  of  scaling 
from  tunnel  to  flight  on  the  rms  pressure. 


SECTION  8 
CONCLUDING  REMARKS 

First  generation  methods  for  predicting  buffeting  dynamic  pressures  and 
the  resulting  structural  dynamic  response  have  been  developed  that  are  simple 
and  usable.  The  fluctuating  pressure  prediction  method  provides  a  mathemati¬ 
cal  model  for  the  overall  rras  levels  and  spectral  shapes  of  the  buffeting 
pressures.  It  is  assumed  that  the  entire  upper  wing  surface  is  affected  by 
fluctuating  pressures  ,and  curves  are  given  relating  the  mathematical  model 
parameters  to  flight  conditions.  The  response  computation  technique  was  made 
appropriate  to  the  uncertainties  in  the  descriptions  of  the  structure  and  the 
buffeting  pressures.  The  assumptions  reduce  the  amount  of  input  data  and 
computational  time. 

Validation  of  wind  tunnel  pressure  measurements  made  on  the  model  was 
accomplished  by  comparison  with  flight  data.  Although  limited  by  the  small 
number  of  flight  transducers,  these  comparisons  showed  good  agreement  in  the 
spectral  shapes.  In  general,  the  model  rms  levels  were  somewhat  higher.  Air¬ 
craft  response  predictions  are  therefore  expected  to  be  somewhat  conservative. 

Comparison  of  predicted  and  measured  dynamic  response  of  the  F-U  airplane 
in  moderate  and  heavy  buffet  showed  reasonable  agreement  in  the  spectral 
shapes  and,  for  the  more  severe  buffet  condition,  in  the  overall  rms  level. 

For  the  cases  of  moderate  buffet,  the  predicted  rms  level  was  higher  then 
that  measured. 

Several  pertinent  conclusions  drawn  from  oho  test  data  include  the 
following; 

o  At  high  angles  of  attack  the  flow  pattern  over  a  highly  swept  wing 
with  a  leading  edge  snag  is  quite  complex.  The  flow  close  to  the 
upper  surface  is  strongly  affected  by  a  vortex  system  consisting  of 
leading  edge,  tip,  and  snag  vortices  causing  abrupt  changes  in  flow 
direction  over  the  surface.  The  regions  of  high  fluctuating  pressures 
seem  to  be  associated  with  these  vertices,  particularly  the  snag  and 
tip  vortices. 

o  K©  discernible  single  "convection”  mechanism  for  the  transport  of  the 
fluctuating  pressures  vm  evident  such  as  the  downstream  convection 
reported  in  studies  of  slender  spacecraft  launch  vehicles.  Distur¬ 
bances  seemed  to  emanate  from  multiple  sources  simultaneously  und 
propagate  in  a  complex  manner. 


o  The  fluctuating  pressure  spectrums  frequently  exhibited  peaks  at  fre¬ 
quencies  believed  to  be  associated  with  vortices.  The  frequency 
associated  with  the  snag  vortex  was  generally  about  half  that  of  the 
leading  edge  and  tip  vortices.  These  frequencies,  scaled  to  the  flight 
vehicle,  were  at  the  higher  end  of  the  vehicle's  primary  structural 
resonances . 

o  Maximum  buffet  intensities  occurred  at  the  high  subsonic  Mach  numbers 
diminishing  abruptly  to  small  values  at  sonic  and  supersonic  speeds, 
o  Maximum  rms  wing  root  bending  moments  due  to  buffet  were  of  the  order 
of  1%  of  the  corresponding  steady  bending  moments, 
o  Maximum  fluctuating  pressure  coefficients  were  generally  of  the  order 

of  ACPrms  =  0,2 * 

The  response  prediction  method  presented  in  this  report  is  based  on 
techniques  of  observation,  measurement,  and  analysis  that  are  considered 
to  provide  a  firm  basis  for  development  effort.  It  therefore  provides  a 
rational  usable  approach  to  the  problem  of  assessing  the  nature  of  the 
fluctuating  pressure  distribution  and  the  structural  dynamic  response  of  a 
highly  swept  fighter  wing. 


SECTION  9 
RECOMMENDATIONS 

This  report  presents  a  first  generation  method  for  predicting  the  dynamic 
response  of  an  airplane  wing  to  input  buffet  pressures.  Simplifying  assump¬ 
tions  have  been  made  in  regard  to  flow  characteristics,  mathematical  modeling, 
and  computational  procedures.  The  simplifications  employed  were  based  upon 
existing  technology,  experience  gained  from  earlier  buffet  research,  and  the 
realization  that  high  levels  of  sophistication  may  not  materially  improve  the 
predictions.  The  method  developed  here  was  used  to  predict  the  response  of 
an  P-4  wing  to  buffeting  pressures,  and  Section  7  presents  the  results  in 
comparison  with  flight  data. 

Additional  refinement  to  the  method  should  be  based  upon  carefully  con¬ 
ducted  evaluation  studies  of  the  possible  areas  of  refinement.  Thorough 
study  will  help  to  insure  that  modifications  do  provide  improvement  and  de¬ 
finite  advantage  over  the  existing  method.  Two  possible  levels  of  refinement 
to  the  prediction  method  presented  herein  are  immediately  available.  The 
first  consists  of  refining  the  existing  mathematical  model  so  that  it  can  be 
readily  extended  to  other  aircraft.  The  second  involves  an  increase  in  the 
complexity  of  the  mathematical  model  based  on  a  better  understanding  of  the 
flow  phenomena.  These  two  levels  of  refinement  are  discussed  below; 

(1)  Extension  to  other  aircraft  of  the  present  buffet  loads  prediction 
method  could  be  enhanced  by  further  wind  tunnel  tests  to  determine 
the  effect  of  variations  in  wing  geometry  on  ouch  mathematical 
parameters  as: 

©  Reduced  frequency,  kt> ,  characteristic  of  the  plan form  vortices 
o  Spectral  amplification  factor,  n ,  corresponding  to  the  char¬ 
acteristic  reduced  frequency 
o  Reduced  spatial  decay  parameter,  6e 
The  variations  in  wing  geometry  might  include: 

©  Taper 
o  Sweep 
©  Camber 
o  Airfoil  shape 

(2)  If  a  comparison  of  measured  data  with  results  of  1  he  application 
of  the  prediction  method  as  described  in  Item  (l)  indicate  that 
more  precise  predictions  are  needed,  then  further  research  might 


be  pursued  as  follows: 

a.  Define  a  systematic  analytical /experimental  attack  on  the 
problem,  building  successively  on  bases  established  in  progres¬ 
sing  from  the  simpler  to  the  more  complex  problem  parameters. 

:  It  is  suggested  that  an  initial  base  of  understanding  might 

be  developed  for  subsonic  buffet  prior  to  progressing  to  tran¬ 
sonic  buffet  along  with  a  progression  from  simple  unswept 
planforms  (without  leading  edge  vortices)  leading  to  the  highly 
swept  and  delta  planforms.  Snag  and  other  planform  anomalies 
should  be  deferred  to  the  later  steps  in  the  development  to 
avoid  obscuring  basic  geometry  effects. 

b.  Analytical  sophistication  should  be  the  minimal  necessary  for 
adequate  representation  at  each  stage  and  increased  only  suf¬ 
ficiently  to  account  for  the  next  level  of  development. 

c .  The  experimental  program  should  parallel  the  analytic  develop¬ 
ment  for  logical  evolution  of  the  mathematical  >pdel. 

Several  specific  topics  discussed  in  thih  report  might  be  included  in 
the  investigation  described  In  (?)  above  or  might  merit  furtner  examination 
independent  of  the  prediction  method.  Tho^r  are: 

o  Identification  of  suitable  tunnels  with  low  turbulence  background 
o  toady  state  flow  conditions  over  thee  wing  'tV.g. »  locations  and  sixes- 
of  vortices,  no  pu  rat  ion  Pubb  ten ,  local  convert  loti  direction;.! ,  regions, 
depurated  flows)  defined  as  fuii'-i. louct  of  flow  erai  lltioji'  ( fiur-stream 
flow  condition  and  mmlui  orientation)  either  fr^i "analysis’ 'or  experi¬ 
ment  .  Oi  l  flow  umi  tuft  ttu'hui qUviii  need  to  InPcuppl ejnsntofi  'Vi th 
additional  investigation;:; ,  such  an 

*  l'-u j'veyt'.  >>r  pi  UM.-utut t>>  profiles  with  nplghl  a'ncvo  the  wing 
■uirfae*1  at  n  number  of  planform  i.«enti on& 

*  durv-y-;  e.f  flow  d  I  reel,  ton  ppefi  via  probe  mounted  tufts  on  and 
■d*cn<>  <  h*>  vo'fntH-  at  the  ‘rune  {■!  ue'orm  toent  iuus  un*id  In. the  pitot- 
staid. *  survey.!.  to  whether  tuft  it;  steady  dr  buffeting. 
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Further  examination  of  buffet  characteristics  "between  M  =  0.9  and  M  = 
1.0  where  there  is  a  sharp  reduction  in  intensity. 

0  Scaling  of  rms  pressure  coefficient  from  model  to  full  scale. 

0  Effect  of  altitude  on  structural  dynamic  response. 

0  Additional  development  of  mathematical  model  for  low  frequency  slope 
in  pressure  power  spectral  densities.  .. 

0  Convection  patterns  on  three-dimensional  wings. 

°  Effect  of  using  average  mode  approach  in  development  of  generalized 
forces. 

°  Effect  on  response  prediction  of  neglecting  correlation  between  areas. 

0  Effect  of  neglecting  normal  coordinate  cross  power  spectral  den¬ 
sities  on  response  prediction. 

0  Provision  for  greater  flexibility  in  the  definition  of  regions. 

°  Optimization  of  the  sizing  of  regions  and  basic  areas. 

The  severity  of  the  buffet  dynamic  loads  problem  should  be  evaluated  on 
a  broader  basis  than  was  possible  with  the  existing  data.  The  limited  exist¬ 
ing  data  from  this  and  other  recent  programs  suggest  that  rms  wing  root  bend¬ 
ing  moments  due  to  buffet  are  a  small  fraction  of  the  static  bending  moment , 
i.e.,  5-10#  more  or  less.  Similar  data  becoming  available  after  this  program 
should  be  collected  and  evaluated  for  a  broader  assessment  of  the  problem 
severity.  Toward  this  end,  future  flight  and  wind  tunnel  buffet  programs 
should  provide  tape  records  of  the  bending  moments  for  subsequent  processing 
to  separate  the  fluctuating  bending  moments  from  the  static  bending  moments. 
Several  different  spanwise  locations  should  be  included  to  help  identify  the 
more  critical  areas.  Collection  of  those  data  might  also  lead  to  a  more 
direct  technique  for  preliminary  design  prediction  of  buffet  loads. 
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Figure  1  Buffet  Loads  Prediction  Method 
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Fluctuating  Pressure  Oscillograph  Traces  for  Configuration  1,  Mach  0.7,  Wing  Angle 


Torsion,  and  Acceleration  Oscillograph  Traces  for  Configuration  1,  Mach  0.7,  Wing  Angle 


Fluctuating  Pressure  Oscillograph  Traces  for  Configuration  1,  Mach  0.9,  Wing  Angle  of  Attack  12° 


Figure  It  Bending,  Torsion,  and  Acceleration  Oscillograph  Traces  for  Configuration  1,  Mash  0.9,  Wing  Angle 


Honeywell  Model  9300B  Automatic  Spectrum  Analyzer 


Diagram  of  Spectral  Data  Reduction  at  MCAIR 
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Figure  15  Root  Mean  Square  Pressure  Coefficient  vs  Percent  Chord  (Configuration  1 

Mach  0.7,  Wing  Angie  of  Attack  12°) 
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Figure  18  Root  Mean  Square  Pressure  Coefficient  vs  Percent  Chord 
(Configuration  1,  Mach  0.9,  Wing  Angle  of  Attack  12°) 
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Figure  20  Root  Mean  Square  Pressure  Coefficient  vs  Percent  Chord 
( Configuration  1,  Mach  1.0,  Wing  Angle  of  Attack  0°) 
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Figure  21  Root  Mean  Square  Pressure  Coefficient  vs  Percent  Chord 
<  Configuration  1,  Mach  1.0,  Wing  Angie  of  Attack  12° ) 
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Figure  23  Root  Mean  Square  Pressure  Coefficient  vs  Percent  Chord. 
{Configuration  5.  Mach  0.9,  Wing  Angle  of  Attack  11°) 
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Figure  24  Root  Meat)  Square  Pressure  Coefficient  vs  Percent  Chord 
(Configuration  8,  Mach  0.9,  Wing  Angle  of  Attack  14°) 


Root  Mean  Square  Pressure  Coefficient  (AC 


.  i 


Steady  B< 


Figure  29  Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.7, 
Wing  Angle  of  Attack  -1  Deg..  Transducer  10) 


Fluctuating  Pressure  Power  Spectra!  Density  (Configuration  1,  Mach  0.7, 
Wing  Angle  of  Attack  -1  Deg.,  Transducer  13) 
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Figure  36  Fluctuating  Pressure  Power  Spectral 
Wing  Angle  of  Attack  -1  De 


Figure  37  Fluctuating  Pressure  Power  Spectral  Density  (Configuration 
Wing  Angle  of  Attack  0  Deg.,  Transducer  13) 


Figure  39  Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.7, 
Wing  Angfe  of  Attack  12  Deg.,  Transducer  1) 


Figure  40  Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.7 
Wing  Angie  of  Attack  12  Deg.,  Transducer  2) 


Fluctuating  Pressure  Power  Spectral  Density  (Configuration 
Wing  Angle  of  Attack  12  Deg.,  Transducer  3) 
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Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.7 , 
Wing  Angle  of  Attack  12  Deg.,  Transducer  11) 


Fluctuating  Pressure  Power  Spectral  Density  (Configuration 
Wing  Angle  of  Attack  12  Deg.,  Transducer  12) 
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Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1 ,  Mach  0. 
Wing  Angle  of  Attack  12  Deg.,  Transducer  13) 
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Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.7, 
Wing  Angle  of  Attack  12  Deg.,  Transducer  16) 


Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.7, 
Wing  Angle  of  Attack  12  Deg.,  Transducer  IS) 
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Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.7, 
Wing  Angle  of  Attack  12  Deg.,  Transducer  21) 
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Figure  61  Fluctuating  Pressure  Power  Spectral  Density  (Conf  iguration  1,  Mach  0.7, 
Wing  Angle  of  Attack  12  Deg.,  Transducer  24) 


Spectral  Density  (Configuration  1,  Mach  0.9, 
:k  8  Deg.,  Transducer  10} 
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Hyure  65  Wing  Locations  at  which  Peaks  are  Evident  in  Pressure  Power  Spectral  Density  Functions 
(Configuration  1,  Mach  0.9,  Wing  Angle  of  Attack  8°) 
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Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.9, 
Wing  Angie  of  Attack  12  Deg.,  Transducer  10) 


Figure  67  Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.9, 
Wing  Angle  of  Attack  12  Deg.,  Transducer  11) 


Figure  68  Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.9, 
Wing  Angle  of  Attack  12  Deg.,  Transducer  13) 
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Figure  70  Wing  l  ocation!  at  which  Peaks  are  Evident  in  Pressure  Power  Spectral  Density  Functions 
(Configuration  1,  Mach  0.9,  Wing  Angle  of  Attack  12°) 
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Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.9, 
Wing  Angle  of  Attack  14  Deg.,  Transducer  13) 


Figure  76  Fluctuating  Pressure  Power  Spectral  Density  (Configuration 
Wing  Angle  of  Attack  12  Deg.,  Transducer  11) 


Figure  77  Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1.  Mach  1.0, 
Wing  Angle  of  Attack  12  Deg.,  Transducer  13) 


Figure  79  Fluctuating  Pressure  Power  Spectral  Density  (Configuration  5,  Mach  0.7, 
Wing  Angie  of  Attack  1 1  Deg.,  Transducer  10} 


Figure  80  Fluctuating  Pressure  Power  Spectral  Density  (Configuration  5,  Mach  0.7, 
Wing  Angie  of  Attack  11  Deg.,  Transducer  11} 


Fluctuating  Pressure  Power  Spectral  Density  (Configuration  5,  Mach  0.7, 
Wing  Angle  of  Attack  11  Deg.,  Transducer  14) 


Fluctuating  Pressure  Power  Spectral  Density  (Configuration  5,  Mach  0.9, 
Wing  Angle  of  Attack  1 1  Deg.,  Transducer  11) 


ower  Spectral  Density  (Configuration  5,  Mach  0.9, 
Attack  11  Deg.,  Transducer  13) 


e  86  Fluctuating  Pressure  Power  Spectral  Density  (Configuration  5,  Mach  0.9, 
Wing  Angle  of  Attack  1 1  Deg.,  Transducer  14) 
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Figure  88  Fluctuating  Pressure  Power  Spectral  Density  (Configuration  5,  Mach  0.9, 
Wing  Angle  of  Attack  14  Deg.,  Transducer  10) 


Figure  89  Fluctuating  Pressure  Power  Spectral  Density  {Configuration  5,  Mach  0.9, 
Wing  Angle  of  Attack  14  Deg.,  T ransducer  1 1 ) 
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Fluctuating  Pressure  Power  Spectral  Density  (Configuration  6,  Mach  0.7, 
Wing  Angle  of  Attack  11  Deg.,  Transducer  10) 


Figure  94  Fluctuating  Pressure  Power  Spectral  Density  (Configuration  6,  Mach  0.7, 
Wing  Angle  of  Attack  1 1  Deg..  Transducer  1 1 1 


Fluctuating  Pressure  Power  Spectral  Density  (Configuration  6,  Mach  0.7, 
Wing  Angle  of  Attack  11  Deg.,  Transducer  13) 


Percent  Chord  (x/c  x  100) 


Figure  98  Ratio  of  Peak  to  Average  Low  Frequency  Power  Spectral  Density  Level  vs  Percent  Chord 
(Configuration  1,  Mach  0.9,  Wing  Angle  of  Attack  8°) 
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Figure  99  Ratio  of  Peak  to  Average  Low  Frequency  Power  Spectr 
(Configuration  1,  Mach  0.9,  Wing  Angle  of  A 


Ratio  of  Peak  to  Average  Low  Frequency  Power  Spectral  Density  Level  (ij) 


Fi«uru  102  Ratio  of  Peak  to  Average  Low  Frequency  Power  Spectral  Density  Level  vs  Percent  Chord 
(Configuration  5.  Mach  0.9,  W i..u  Angle  of  Attack  11°) 


Figure  106  Fluctuating  Pressure  Normalized  Cross  Power  Spectral  Density  (Phase) 
(Configuration  1,  Mach  0.7,  Wing  Angle  of  Attack  12  Deg.,  Transducers  6  and  10) 


gure  W?  Fluctuating  Pressure  formalized  Cross  Power  Spectral  Density  (Magnitude) 
(Configuration  1,  Mach  0.7.  Wing  Angle  of  Attack  12  Deg.,  Transducers  12  and  18) 
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:igure  110  Fluctuating  Pressure  Normalized  Cross  Power  Spectral  Density  (Phase) 
(Configuration  1,  Mach  0.9,  Wing  Angle  of  Attack  8  Deg.,  Transducers  6  and  10) 


Figure  114  Fluctuating  Pressure  Normalized  Cross  Power  Spectral  Density  (Phase) 
(Configuration  1,  Mach  0.9,  Wing  Angle  of  Attack  12  Deg.,  Transducers  6  and  10) 
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gure  1 19  Fluctuating  Pressure  Normalized  Cross  Power  Spectra!  Density  (Magnitude) 
(Configuration  1.  Mach  0.9,  Wing  Angle  of  Attack  14  Deg.,  Transducers  12  and  18) 
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Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.9, 
Wing  Angle  of  Attack  12  Deg.,  Transducer  27) 
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Fluctuating  Pressure  Power  Spectral  Density  (Configuration  1,  Mach  0.9, 
Wing  Angle  of  Attack  8  Deg.,  T ransducer  21 ) 
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gure  131  Fluctuating  Pressure  Normalized  Cross  Power  Spectral  Density  (Magnitude) 
(Configuration  1,  Mach  0.9,  Wing  Angle  of  Attack  12  Deg.,  Transducers  21  and  27) 
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gure  135  Comparison  of  Model  and  Full  Scale  rins  Pressure  Coefficient  Variation  with 
Angle  of  Attack  at  Mach  0.8, 90  Percent  Chord.  33  Percent  Semispan 
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Figure  136  Comparison  of  Model  and  Full  Scale  mis  Pressuiu  Coefficient  Variation  with 
Angle  of  Attack  at  Mach  0,8,  95  Percent  Chord,  33  Percent  Semispan 
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Figure  137  Comparison  of  Model  and  Full  Scale  rms  Pressure  Coefficient  Variation  with 
Angle  of  Attack  at  Mach  0.9, 90  Percent  Chord,  86  Percent  Semispan 


- T - r - - r - 


Configuration  t 


O  f  till  Huaia.  Altitude  » 2ti,0tMi  ft.  q  •  4WJ  pst 


A  Full  ScaW,  Alliluda  -  dli.OOO  ft,  q  280  p»i 


Q  Mqtlal,  q  •<  640  Pit 


Mapraiantdlivti  Cuiva  Itn  all  Frill  Scott-  t>«ia 


w 

til 


VVhhj  A««(Ii>  of  Attack  lit  i,  tlitg 


Figure  1 38  Comparison  of  Model  and  Full  Scale  rms  Pressure  Coefficient  Variation  with 
Angle  oi  Attack  at  Mach  0.9, 90  Percent  Chord,  78  Percent  Semispan 
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Ratio  of  Peak  to  Average  Low  Frequency  Power  Spectral  Density  Level  (i?) 

Figure  150  Buffet  Pressure  Prediction  Parameter  X  as  a  Function  of  Parameter 


Comparison  of  Theoretical  and  Measured  Power  Spectral  Density  Spectra 
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Figure  162  Nondimensionul  Frequencies  Associated  with  Leading  Edge  Tip, 
and  Snag  Vortices  (Continued) 
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Figure  152  Nondimensional  Frequencies  Associated  with  Leading  Edge,  Tip, 
and  Snag  Vortices  (Concluded) 
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Figure  1b3  Model  for  Rotio  of  Peak  to  Average  Low  Frequency 
Power  Spectral  Density  Level  (»)  (Continued) 
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Fiyuro  155  Example  of  Division  of  Wing  Plunform  into  Regions  and  Areas 


Acceleration  Power  Spectral  Density,  g^/Kz 
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Figure  159  Effect  of  Damping  Ratio  on  Predicted  Acceleration  Response 


Acceleration  Power  Spectral  Density,  g  /Hz 


VISUALIZATION  STUDIES 

Methods  of  Analysis 

Information  available  from  the  wind  tunnel  test  for  use  in  flow  visuali¬ 
sation  studies  included  l6mm  color  movies  and  black  and  white  still  photos 
of  the  oil  and  tufts  on  the  top  surface  of  the  wing. 

Color  movies  of  the  wind  tunnel  tests  were  valuable  tools  in  the  ana¬ 
lysis.  The  films  were  run  at  slow  speeds  and  could  be  stopped  for  more 
detailed  study  of  any  given  frame.  From  tne  flow  visualization  movies, 
flow  charts  were  sketched  and  observations  were  recorded. 

The  tufts  were  classified  by  degree  of  oscillation  and  flow  direction 
r.s  follows: 

1)  STABLE  -  clearly  definable  flow  direction,  little  or  no  oscillation  - 
indicates  undisturbed  airflow. 

2)  MILD  -  clearly  definable  flow  direction,  base  of  tuft  stable,  end 
of  tuft  in  a  bounded  oscillation. 

3)  MODERATE  -  definable  flow  direction;  entire  tuft  oscillates  in 
a  bounded  area. 

h)  RANDOM  -  indefinable  flow  direction;  tuft  oscillates  with  no 
definite  pattern  -  indicates  separated  airflow. 

The  tuft  flow  directions  were  marked  on  scale  drawings  of  the  wings  as 
were  lines  separating  different  flow  directions. 

fhv  <41  vug  analyzed  in  a  (similar  manner.  ‘.Sketches  were  made  from 
the  movies  shaving  flow  directions,  pattern,  and  accumulations  of  the  oil. 

Mm-  tu  flow  directed  aft  and  outboard  was  classified  as  attached.  Areas 
wiv-rc  Ur-  oil  collected  were  classified  as  separated  regions.  Areas  where 
the  oil  has  h  forward  flow  direction  were  also  classified  as  separated. 

otiU  photos  of  selected  conditions  were  used  for  comparative  a*.uiysis. 
Flow  patterns  and  general  trends  were  correlated  between  oil  and  tufts. 

•Vitoria  used  in  analyzing  movie:*,  were  also  used  for  still  photos.  The 
Vr.rev  f  oscillation  of  the  tufts  was  determined  by  the  size  of  the  blur 
■>:'  each  tuft,  and  the  flow  direction.  It  appears  that  as  the  tufts  oscilluted, 
‘heir  bujk  tuid  forth  motion  against  tin*  surface  of  the  wing  polished  the 
air.t  and  formevl  an  arc-shaped  area  which  reflected  light  into  the  cameras. 

A*,  first  glance,  these  light  areas  may  be  thought  to  be  the  tufts  in  oscil- 
.«*  1  .  However,  tufts  in  oscillation  appear  as  a  uniform  blur,  not  as 

w.  'u  vi'-t,  the  tuft  clearly  visible. 
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In  several  cases,  the  tufts  and  the  oil  appeared  to  be  indicative  of 
different  types  of  flow.  Differences  in  the  tuft  and  oil  patterns  may  be 
due  to  different  streamlines  at  and  near  the  surface  of  the  wing.  The 
tufts  are  influenced  by  a  layer  as  thick  as  the  tufts,  containing  several 
streamlines.  This  gives  the  tufts  three  dimensions  in  which  to  react. 

The  oil  reacts  only  two  dimensionally  because  it  is  affected  only  by  the 
air  on  the  surface  of  the  wing.  In  some  cases  with  leading  edge  separation, 
the  oil  flow  flowed  forward  from  the  oil  ports  and  did  not  reach  the  area 
aft  of  the  oil  ports.  In  such  circumstances ,  the  tufts  were  the  better 
visual  metnod  of  determining  the  type  of  airflow  on  the  aft  portion  of 
the  wing. 

On  the  F-4  wing,  there  is  a  system  of  leading  edge  vortices  which  are 
responsible  for  many  of  the  tuft  and  oil  patterns.  Their  boundaries  can 
be  approximated  by  the  changes  of  directions  of  the  tufts,  and  by  the  cir¬ 
culatory  behavior  of  the  oil.  The  F-4  wing  is  like  f  o  wings  in  one,  each 
with  a  leading  edge  vortex.  See  Figure  A-l.  The  outer  wing  Juts  for¬ 
ward  of  the  inner  wing,  creating  a  snag  which  generates  another  vortex 
aft  along  the  fold  line.  Another  vortex  (the  tip  vortex)  comes  from 
beneath  the  wing  and  rolls  over  the  wing  tip.  This  complex  system  of 
vortices  is  a  partial  cause  of  leading  edge  separation,  regions  of  oil 
circulation,  different  flow  directions  of  the  tufts,  and  the  accumulation  of 
oil  in  regions  of  low  pressure. 

Results 

Variation  of  Flow  Characteristics  with  Angle  of  Attack  -  These  studies  wer 
run  for  Configuration  1  at  Mach  0.9  and  Reynolds  Humber  of  1.2b  x  10  /ft. 

At  aw  -  i»°,  see  Figure  A-2 ,  the  oil  flows  smoothly  aft  over  the  largest 
portion  of  the  wing.  There  is  some  collection  of  the  oil  along  the  leading 
edge  on  the  outboard  panel,  indicating  separation  which  is  probably  caused 
by  the  leading  edge  vortices.  The  tufts  on  the  area  corresponding  to  the 
oil  collection  are  directed  outboard  and  are  parullel  to  the  loading  edge. 

The  tu,.ts  on  the  outboard  panel  ure  in  mild  to  moderate  oscillation.  The 
'ufts  on  the  inner  wing  range  from  stable  to  mild  oscillation. 

At  «w  »  8°  (Figure  A-3),  the  oil  over  most  of  the  wing  flows  outboard 
arid  aft,  indicating  smooth  attached  airflow.  The  leading  edge  vortices 
cause  the  oil  to  flow  outboard  along  the  leading  edge  and  collect  in  the 
areus  of  separated  flow.  The  abrupt  changes  in  the  flow  directions  of  the 


tufts  are  evidence  of  the  leading  edge  vortices.  Tufts  on  the  outboard 
panel  approach  random  motion,  indicating  separation.  The  only  region  of 
stable  tufts  lies  next  to  the  fuselage  close  to  the  >•*,  edg- implying 

undisturbed  airflow. 

At  av  =  12°  (Figure  A-t),  the  oil  is  swept  outboard,  parallel  to  the 
leading  edge,  until  it  hits  the  fold  line  and  is  deflected  slightly  aft. 

The  tufts  show  the  effect  of  the  snag.  The  vortex  generated  by  the  snag 
deflects  the  innermost  leading  edge  vortex  and  sweeps  off  the  outboard 
panel  (Figure  A-l).  The  random  motions  of  the  tufts  on  the  ou+er  two  thirds 
of  the  wing  indicate  that  the  airflow  is  detached. 

'  -  =  ih°  (Figure  A-5) , the  oil  collects  on  the  inboard  leading  edge, 

t-  ' 

■. 1  /ug  1  f'  .  and  along  a  line  aft  of  the  outboard  leading  edge 

lap,  showing  the  vo»-t. !-*►•*  and  indicating  separation  in  the  respective  areas. 
Tufts  in  random  motion  indicate  separation  over  a  large  percentage  of  the 
wing.  Tufts  along  the  leading  edge  show  that  the  flow  is  directed  inboard. 

This  is  caused  by  the  vortices  and  the  fact  that  at  high  angles  of  attack 
the  pressure  along  the  leading  edge  Is  lower  inboard.  On  the  tufted  wing, 
the  forward  boundary  of  the  leading  edge  vortex  is  indicated  by  the  line 
of  demarcation  (Figure  A-5)  extending  aft  and  outboard  from  the  inboard 
leading  edge.  Inboard  of  the  line,  the  tufts  point  outboard,  but  outboard 
of  the  line  they  are  directed  inboard. 

At  aw  =  20°  (Figure  A-6),  the  oil  shows  that  on  the  forward  inboard  por¬ 
ts  on  of  the  wing  the  flow  moves  forward  and  inboard,  aft,  and  then  outboard. 

This  on  is  caused  by  the  innermost  leading  edge  vortex  and  the  pressure 
uifitr  i.  but  ion  along  the  leading  edge.  The  center  of  circulation  can  be 
approximated  from  the  flow  directions  of  the  tufts.  On  the  inboard  leading 
edge  a  spinning  tuft  is  found  which  corresponds  with  the  area  of  circulating 
'■.Ll.  The  tufts  indicate  separation  over  the  entire  wing.  At  this  high 
iiigle  of  attack,  the  leading  edge  vortex  i.-  above  the  surface  of  the  wing; 

:  uwe/er,  the  tufts  are  still,  affected  and  indicate  the  location  of  its 
.  .••ward  boundary. 

As  i  i. .  wing  nng.l  o  of  (■<  )  i  acre  as  ess ,  the  percentage  of  turbulent 

W 

a."..  .  ir.  si  ••ated  :.y  the  tufts,  increases  and  the  leading  edge  vortices  move 
;;  i  -:rd.  Areas  of  separation  are  relatively  stationary  in  their  location, 

•■■is  mere  ..Sr?  in  size  and  ir.t:  nslty  as  u  gets  larger. 
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Variation  of  Flow  Characteristics  with  Mach  "umber  -  These  studies  were 
:  an  for  Configuration  1  at  =  12°  and  Reynolds  Number  of  1.25  x  10*7  ft. 

At  M  =  0.7  (Figure  A-7),  the  oil  on  the  inner  wing,  m-.-ven  mu  and  out¬ 
board,  demonstrating  smooth  attached  flow.  The  collection  of  oil  along  the 
leading  edge  indicates  separation.  The  oil  on  the  remainder  of  the  outboard 
panel  is  swept  off  by  the  system  of  vortices.  The  random  motions  of  the  tufts 
on  the  outboard  panel  show  this  entire  area  is  separated.  The  tufts  on 
the  remainder  of  the  wing  are  in  mild  to  moderate  oscillation. 

At  M  ■  0.9  (Figure  oil  no  longer  c< luects  on  the  leading  edge. 

The  tufts  indicate  that  tue  area  of  separated  flow  nas  increased. 

At  M  a  1.0  (Figure  A-8),  the  oil  begins  to  circulate  in  two  areas: 

l)  on  the  inboard  panel  behind  the  leading  edge  flap,  and  2)  around  the  fold 

line,  encompassing  the  entire  outboard  panel.  These  oil  patterns  ar.  caused 
by  a  combination  of  the  leading  edg*  vortices,  the  snag  vortex,  separation, 
a?  a  the  t  •~uro  distribution  over  the  wing.  The  separated  area  shown  by 
the  tufts  has  decreased  from  that  at,  a  -  omull  perc-'-’tag**.  This 

trend  seems  to  indicate  a  critical  point  between  a. 9  and  1.0  (0.9  <  M  ’  <  1.0). 

»  ““  viil  X  "~ri 

At  M  =  1.2  (Figure  A-9),  the  oil  present  on  the  center  of  th°  wing 

collects  on  the  leading  edge,  along  a  possible  veuk  shock.  The  •  "emula¬ 

tion  of  oil  on  the  wing  tip  is  due  to  the  leading  edge  vortices,  including 
the  snag  vortex  and  the  tip  vortex.  The  oil  on  the  outboard  panel  gathers 
along  a  shock  line  parallel  to  the  trailing  edge.  Reparation  indicated 
by  the  tufts  is  limited  to  the  periphery  of  the  outer  wing  panel.  The 
tufts  in  random  motion  at  the  trailing  edge  correspond  with  the  area  aft 


ot  the  shock  line  in  the  oil.  The  boundaries  of  the  leading  edge  vortices 

are  definable  by  the  changes  in  the  flow  directions  of  the  tufts. 

At  a  constant  angle  of  attack,  the  area  of  separation  increases  as 

M  increases  to  a  point  between  0.9and  1.0,  M  .  Above  M_ the  area  of 

Uni L  CnIT 

o^puration  decreases,  at  least  no  to  M  =  1.2. 

/aviation  of  Flow  Characteristics  witm  reyrmlds  number  —  These  studies 
wore  run  for  Configuration  1  at  Mach  0.9. 

At  «v  -8°,  Re  =  1.25  x  10  /ft.  (Figure  A- 3)  ,  the  leading  edge  vortices 
cause  the  oil  to  collect  along  the  leading  edge  on  the  outboard  panel , 
ami  just  inboard  of  the  snag.  The  tufts  on  the  outboard  panel  are  in 
moderate  oscillation.  The  area  of  stable  tufts  is  small  and  is  inboard 
ah  :.g  the  leading  edge.  This  type  of  flow  is  typical  for  a  highly  swept 


these  conditions. 
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At  a  8°,  Re  3.75  x  10  /ft.  (Figure  A-10)  ,  there  is  no  accumulation 
of  oil  on  the  inboard  panel,  and  the  oil  on  the  outboard  panel  is  closer 
to  the  stag.  This  might  imply  ‘•hat  an  '**  j*v:r“-  *»•*  t  the  vor* ices  -ppoar  ic 
shift  slightly  inboard.  Tiv  cine  >,ird  panel  in  similar  to  that  with  the 
lower  Reynolds  number,  but  the  stable  area  of  tufts  has  increased.  This 
stabilizing  effect  at  the  higher  Re  may  be  due  to  the  thinning  out  of  the 
viscous  layer. 

I't  aw  =  12°.  a  comparison  of  Figures  A-)*  and  A-ll  shows  no  noticeable 
Reynolds  number  effect  in  either  the  tufts  or  the  oil.  An  examination  of 

>thcr  photos  has  shown  no  noticeable  Reynolds  number  effect  for  a  from 

w 


•r,i  1  ,  ■  f  ••  studies  suggest  that  increasing  Reynolds  number  tend 

increase  the  ;  :e  somewhat  at  the  V.wer  angles  of  attack 

near  buffet  onset  but  hue  :m  noticeable  effect  a*  angles  of  attack  well 
beyond  onset .. 


INSTRUMENTATION  STUDIES 


General 

The  instrumentation  data  used  here  to  analyze  the  flow  characteristics 
consisted  of 

o  Aerodynamic  coefficient  data,  C ^ 
o  rms  wing  root  bending  moment, 

o  Steady  pressure  coefficients  from  static  pressure  taps,  C 

* 

o  rms  pressure  coefficients  from  high  frequency  pressure  transducers, 

^Prms  ’ 

These  studies  probe  into  the  flow  characteristics  at  M  =0.7  and  M  =0.9 

and  are  divided  generally  into  two  categories:  buffet  onset  and  post-onset. 

Buffet  Onset. 

Measured  lift  coefficient,  CT ,  versus  angle  of  attack  is  plotted  in. 
Figure  A-12  for  M  =0.7  and  Figure  A-13  for  M  =0.9.  Buffet  onset,  a^,  is 
identified  on  these  curves  by  the  departure  from  linearity,  i.e.  =  9° 
at  M  =  0.7  and  on.,.  =  7°  at  M  =0.9.  The  exert  point  of  departure  from  linea- 
rity  may  be  subject  to  some  uncertainty  because  the  departure  is  not  abrupt, 
as  in  Figure  A-12,  This  is  frequently  the  case,  and  the  location  of  buffet 
onset— to  a  degree— depends  on  the  analyst's  interpretation  of  the  curve 
connecting  the  test  points. 

Another  indicacor  of  buffet  onset  is  the  wing  root  rms  bending  moment. 
This  data  is  presented  in  Figures  A-lU  and  A-15  for  M  =  0.7  and  0.9 ,  respec¬ 
tively.  Here,  buffet  onset  is  identified  with  a  sudden,  sha;  .  rise  in  rms 
from  a  basically  constant  level  below  onset.  This  indicator  presents  a 
rather  good  identification  of  buffet  onset  ,an_  =  8°  at  M  =0.7  and  at  =  7° 
at.  M  =0.9.  It  is  less  dependent  on  the  analyst's  interpretation  than  the 
onset  identification  based  on  CT  . 

XJ 

Static  pressure  (C^)  divergence  is  often  .1  as  an  indicator  of 

buffet  onset.  Such  data  is  presented  in  Figures  A- lb  and  A- 17  for  t>.  =0.7 

and  0.9  respectively  for  the  trailing  edge  {'■)'■  t%  chord)  static  pressure  taps. 

In  Figure  A- 2  0,  for  M  =0.7,  the  first  indication  of  flow  separation  is 

at  a  =  5°  where  the  tip  trailing  edge  static  pressure  shows  a  sudden 
w 

.crease  in  C  from  its  nominal  value  at  the  lower  angles.  At  this  angle 

of  at. tack,  a  =  5°,  the  static  pressures  at  the  other  trailing  edge  span- 


-,o  locations  indicate  that  flow  separation  has  not  occurred  at  these 

however,  at  a  =  8°  all  the  other  trailing  edge  static  pressures 
vv 
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diverge  simultaneously.  From  this  information  it  is  inferred  that 

o  A  local  area  of  flow  separation  initially  occurs  in  the  vicinity 
of  the  tip  trailing  edge  at  =  5° 
o  When  a  reaches  8°  i separation  suddenly  spreads  to  encompass 

V 

the  whole  trailing  edge. 

For  M  =  0.9, Figure  A-17  shows  a  similar  pattern  except  that  flow  separation  t 
over  the  whole  trailing  edge  now  occurs  at  ay  *  6°  instead  of  8°. 

The  validity  in  using  static  pressure  divergence  as  an  indicator  of 
local  flow  separation  (or  at  least  local  pressure  fluctuations)  is  support'd 
by  the  rms  fluctuating  pressure  data,  ACpj^g,  shown  in  Figures  A-l6  and 
A- i 9  V  ~  o/fand  0.9 respectively.  These  show  a  ri3e  in  fluct  ;*ting 
j':ff  ssure  at  the  tip  trailing  edge  at  o  *  5°  for  M  *  O.Tand  o  *  1»°  for  • 
j-i  *  o,9  reasonably  corresponding  to  the  static  pressure  divergences  for 
the  same  location.  They  further  show  fluctuating  pressure  rises  for  the 
other  trailing  edge  locations  at  about  ay  *  7°  for  H  ■  o.7and  »  6°  for 
M  =  o.9  again  reasonably  correr ponding  to  the  static  pressure  divergences 
for  the  same 'locations. 

Buffet  onset  indications  based  on  and  res  wing  root  bending  moment 
are  ird Icative  of  the  gross  flow  picture  over  the  wing  while  the  static 

pr-rsures,  C  ,  are  indicative  of  local  flow  conditions  existing  near  the 

P 

statin  prec jure  taps.  As  a  consequence, static  prt-ture  divergence  is 
Ludicr.r. i-e  of  the  initiation  and  growth  of  local  areas  of  flow  separation, 
w.' :  i  e  Cr  and  rms  bending  moment  divergences  are  indicative  of  sufficient 
V/- in  local  flow  separation  to  be  detectable  on  a  gross  scale.  The 
i r  i.i  what  is  loosely  defined  as  buffet  onset.  Even  though  a  small 
o’  local  flow  separation  (tip  trailing  edge)  is  indicated  at  =  5° 
for  1  :‘.h  M  =0.7  andO.9,  buffet  enset  is  not  considered  to  have  occurred 
until  i.  reaches  8°  or  9°  at  M  =0.7  and  7°  at  M  =0.9. 

A 

:  •••••••  uswt 

A--,  mentioned.  elsevher*  in  this  r  .port ,  the  F-l  wing  with  its  high 
-.v'.-ep ,  nigh  taper  and  snag  gives  rise  to  ccr.pl ex  flow  patterns  which  are 
•i  :  ‘nn*  at  h '  yh  -.'  .  lee.  i  .oh.  These  studies  attempted  to  probe  into 

L*  •• :  at  •  between  ’he  vortex  systems  and  regions  of  fluctuating 
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Static  pressure  data,  C  ,  is  presented  in  Figures  A-20  and  A-21  for 
M  =  0. land  0.9 respectively.  These  show  spanwise  Cp  variations  for  50$, 

10%,  and  90 %  chord  locations.  The  pressure  peaks  (actually  negative  C 
peaks)  are  believed  to  be  associated  with  vortices  emanating  from  the  tip 
and  from  the  snag.  The  snag  vortex  is  quite  evident  at  the  50$  and  70$ 
chord  locations  while  the  tip  vortex  is  most  evident  at  the  90$  and  70$ 
chord  locations.  ’The  tip  vortex  is  never  evident  forward  of  70$  chord, 
while  the  snag  vortex  effect  tends  to  move  downstream  and  inboard  at  the 
higher  angles  of  attack. 

Rms  fluctuating  pressure  data,  ACp  ,  is  presented  in  Figures  A-22 
and  A-23  for  M  —  0,7&*id  0.9 respectively .  Although  not  absolutely  conclu¬ 
sive,  if  this  data  is  compared  to  the  corresponding  static  pressure  data  of 
Figures  A— 20  and  A-21,  there  is  a  hint  that  the  maximum  fluctuating  pressures 
seem  to  occur  on  either  side  of  the  vortex  cores.  This  suggests  a  model  of 
the  flow  pattern  as  sketched  below. 


Root  Snot)  Tip 


Spanwise  Variation  of  Pressure  Coefficients 


Figure  A  'J  Flow  Visualization  (Oil  anti  Tufts) 

Configuration  1,  Mach  0.9,  Re  1.25  x  106/ft,  Wing  Angle  of  Attack  (uw)  “  8° 
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Figure  A-6  Flow  Visualization  (Oii  and  Tufts) 

Configuration  1,  Mach  0.9,  Re  -  1.25  x  lofyft,  Wing  Angle  of  Attack  (u^,)  -  20° 


}>M 

$ 
+  ' 

' 


251 


k1  WiiiiMii^i  j  *• 


J  ■iinyTim.  . 


1  ‘  1 


Figute  A-10  Flow  Visualization  (Oil  and  Tufts) 

Configuration  1,  Mach  0.9,  Re  “  3.75  x  lofyft,  Wing  Angle  of  Attack  (aw)  =  8° 


Inner  Wing  Bending  Moment  Root  Mean  Square  T!  Inner  Wing  Bending  Moment  Root  Mean  Square 

Fluctuations  (c..  ).  in.-lb  c  Fluctuations  (<7««  ),  in.-lb 
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re  A-14  rms  Inner  Wing  Bending  Moment  vs  Wing  Angie  of  Attack  for 
Configuration  1,  Mach  0.7 


Figure  A-'i5  rms  Inner  Wing  Bending  Moment  vs  Wing  Angle 
of  Attack  for  Configuration  1,  Mach  0.9 
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Figure  A-18  Root  Mean  Square  Pressure  Coefficient  vs  Wing  Angle  of  Attack, 

Configuration  1,  Mach  0.7 
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Figure  A-23  Root  Mean  Square  Pressure  Coefficient  vs  Percent  Semispan. 
Configuration  1,  Mach  0,9 
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